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Foreword 


It is often said that scientific texts quickly become obsolete. 
Why are the Pauli lectures brought to the public today, when 
some of them were given as long as twenty years ago? The 
reason is simple: Pauli’s way of presenting physics is never 
out of date. His famous article on the foundations of quantum 
mechanics appeared in 1933 in the German encyclopedia 
Handbuch der Physik. Twenty-five years later it reappeared 
practically unchanged in a new edition, whereas most other 
contributions to this encyclopedia had to be completely re- 
written. The reason for this remarkable fact lies in Pauli’s 
style, which is commensurate to the greatness of its subject 
in its clarity and impact. Style in scientific writing is a quality 
that today is on the point of vanishing. The pressure of fast 
publication is so great that people rush into print with hur- 
riedly written papers and books that show little concern for 
careful formulation of ideas. Mathematical and instrumental 
techniques have become complicated and difficult; today most 
of the effort of writing and learning is devoted to the acquisi- 
tion of these techniques instead of insight into important 
concepts. Essential ideas of physics are often lost in the dense 
forest of mathematical reasoning. This situation need not be 
so. Pauli’s lectures show how physical ideas can be presented 
clearly and in good mathematical form, without being hidden 
in formalistic expertise. 

Pauli was not an accomplished lecturer in the technical sense 


vil 
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of the word. It was often difficult to follow his courses. But 
when the sequence of his thoughts and the structure of his 
logic become apparent, the attentive follower is left with a new 
and deeper knowledge of essential concepts and with a clearer 
insight into the splendid architecture of reason, which is theo- 
retical physics. The value of the lecture notes is not diminished 
by the fact that they were written not by him but by some of 
his collaborators. They bear the mark of the master in their 
conceptual structure and their mathematical rigidity. Only 
here and there does one miss words and comments of the mas- 
ter. Neither does one notice the passing of time in his lectures, 
with the sole exception of the lectures on field quantization, 
in which some concepts are formulated in a way that may 
appear old-fashioned to some today. But even these lectures 
should be of use to modern students because of their compact- 
ness and their direct approach to the central problems. 

May this volume serve as an example of how the concepts of 
theoretical physics were conceived and taught by one of the 
great men who created them. 


Victor F. Weisskopf 


Cambridge, Massachusetts 


Preface 


When Pauli taught optics to the students of ETH at Ziirich, 
the main interest of this subject, as of mechanics, was to pre- 
pare the students for the ideas of wave mechanics. In this re- 
spect the section in this book on Hamilton’s theory and its 
application to anisotropic media in the section on crystal optics 
are particularly interesting. But otherwise, optics was at Pauli’s 
time a subject of past glory. 

This has radically changed with the advent of the laser, 
which has again put optics on the frontier of science. Subjects 
such as nonlinear optics and holography, which this year won 
a Nobel Prize, have opened up entirely new technologies for 
which a good knowledge of old-fashioned linear classical optics 
is a necessary preliminary. 

This is a fascinating lesson about the unpredictability of the 
evolution of science. And it is perhaps the best justification 
for this English translation. For Pauli was not so much inter- 
ested in detail as in a critical and logical exposition. This 
makes the present lectures a concise and rewarding introduc- 
tion to the subject. 

The German notes prepared by A. Scheidegger and pub- 
lished in 1948 were not without flaws. This was the reason that 
Pauli had encouraged a second version of notes prepared by 
P. Erdés, now Professor at the University of Florida at Tala- 
hassee. If for this translation the first version was chosen after 
ali it was because it reflects better the original spirit and com- 
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pact style of Pauli. The flaws have been eliminated and some 
precision has been added in the comments of the appendix. 
This process of clarification has made the work of the trans- 
lators worth special mention. As editor I am quite pleased to 
see that the hard work of all of us has produced an improved 
version of a course I always liked. 


Charles P. Enz 


Geneva, 2 November 1971 


Pauli Lectures on Physics: 
Optics and the Theory of Electrons 


Chapter 1. Geometrical Optics 


1. FERMAT’S PRINCIPLE 


If there are two homogeneous, isotropic media, M, and M,, 
which are characterized by their indices of refraction, n, 
and ,, then Snell’s law holds: 


Ny sin DO, = Ne sin Ds . pea] 


Figure 1.1 


For n,/n,>1, a8 0, is increased, a point is reached beyond 
which there exists no real solution for ?,. At the critical 


1 
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angle defined by 


- 92) _ ts 

sin #; Bead 
the refracted light propagates parallel to the interface be- 
tween M, and M,. Experiment shows that for 3,>0”, 
total reflection occurs. However, some reflection always 
occurs, whether total or not. In addition to the refracted 
ray, there is also a reflected ray; that is, the original ray 
is split into two. For the reflected ray, we have the law 


O=%,. fie] 


The refracted ray as well as the reflected ray lie in the 
plane of incidence. 

That the index of refraction n is independent of the angle 
of incidence @ is characteristic of the isotropy of the me- 
dium. On the other hand, it should be well noted that 
depends upon the color of the light. In vacuum, by defi- 
nition, n=1 for all colors. Thus, ” is greater than unity 
in material media. 

The law of refraction as well as that of reflection can be 
characterized by means of an extremal property. If we 


Figure 1.2 
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denote the product of index of refraction and light path 
by the optical path JL, 


L=mP,Q+nP.Q, [1.3] 


then the principle formulated by Fermat states that a light 
ray always travels in such a way that 


SL =0. [1.4] 


This law is identical with the laws of refraction and reflec- 
tion given above. 


Proof. a. Law of Refraction 
From Fig. 1.2 we see that 


L = n,a,/cos 3,+ n,a,/c0s 2, , 


A, A, =a, tan#?, + a, tand?, = constant (auxiliary condition) . 


Figure 1.3 


Using the method of Lagrangian multipliers, we obtain 
dL+A18(A,A,) = 0 
for arbitrary variations of #, and #,. Thus, 


e- (* sin 3, al 


nN, Sin P, il 
+4 Tae tae 
cos? 3, cos? 3, 


cos? #, ZA cos? #, 


) 30-0, ( )38.=0, 
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so that 
n, sin 3,+A=0, 


n,8in #,+-A=0. 
Eliminating 4 from these equations yields 
n, sind, = n, sin’, . 


b. Law of Reflection 
If P* is the image of P, then it is seen without com- 
putation (see Fig. 1.3) that 


P.O +O’ P*> 2,0 ae. 


The optical path is, therefore, minimal. 

For curved surfaces, the optical path need not be min- 
imal; it need only be extremal. Later, we shall see that 
since the velocity of light is proportional to 1/n, then the 
optical path is proportional to the transit time of the ray. 
Thus, according to Fermat, this quantity is to be extremal. 
However, the concept of transit time becomes meaningful 
only when considering wave optics. 

If we now consider a medium with a continuously vari- 
able index of refraction n= n(a, y,2), then we can easily 
formulate Fermat’s principle for this case too; correspond- 
ing to Eq. [1.4], we can state the principle as follows: the 
ray travels between the fixed points P, and P, in such a 
manner that 


s5S=ds{nds=0 _ (P,, P, fixed). [1.5] 
If uw is the parameter characterizing a given curve, then 


s=Jo|3(c0) 


t 
du, 
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where 
L; = £;(u) + $2,(u) (i= 1, 2, 3) 


is a neighboring curve, and P, corresponds to a(u;). Then, 
dx, d da,\? 


on fe( Cama) end aw [a(n 


Here, we have used the fact that the differentiations § 
and d are independent, so that Sdv=dédw. If the second 
term in the last formula is integrated by parts, then there 


follows 
dx,\?}-* 
p2 (=) | ) 52,0U 


ss f3 (2 (ce) Tas (mes 
+ Enea) | 


Since by assumption P, and P, are fixed, then 3$a|, 
= $a|,,=0, and the last term vanishes. Then, however, 
58 can vanish for an otherwise arbitrary choice of 5x, only 
if the integrand in the above formula vanishes identically. 
If the parameter wu is chosen equal to the arc length s 
along the extremal curve, and if vector notation is used, 
then this leads to the condition 


Ps 


OL, 


P, 


dds 
grad n = aE (n aa) : [1.6] 


It might be thought that in the whole derivation the 
parameter wu was superfluous in that it could be set equal 
to s right from the start. This, however, cannot be done, 
since while one can indeed choose u= s along the extremal 
curve, this cannot at the same time be done along the 
varied comparison curve. 

If the variation of S corresponding to a variation of the 
endpoints P, and P, is now investigated, subject to the 
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condition that the curve between P, and P, should remain 
an extremum, then there follows 


sS=n (3 “ba 


P, 


[1.7] 


Py 


Here too, Fermat’s principle states only that the optical 
path [nds is to be an extremum; what kind of an ex- 
tremum is not specified. This latter question is of no 
consequence at all in physics. 

If Eq. [1.6] is again decomposed into components, then 
it should be noted that three independent equations do 


not result. Since 
dw,\? 
oa a 
3) a 


a relation must exist between the equations. If one takes 
the scalar product of Eq. [1.6] with da/ds, then there 
follows 


<* gradn = 


ds de de de de’? 


dx d n oe\ _ A dx dn dx d?x 
~ ds ds ie 


so that (dx/ds)-gradn = dn/ds, which is clearly trivial geo- 
metrically. From this it follows that only the vector 
components have physical significance. When expanded, 
Eq. [1.6] yields 


“de ds ° dat” 


One of the well-known Frenet formulas is d?x/ds? = e/R, 
where e is a unit vector in the osculating plane perpen- 
dicular to dx/ds, and £& is the radius of curvature. It then 
follows that 


(grad n), = 


ae 


or [1.8] 


bl® bls 


{grad (log n)}, = 


In the special case where n is a function of only one co- 
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ordinate, for example n(x), then 


d/ dy d/ dz 
as" Ge) =° as") =o 
so that 
dy dz 
n aa constant , artes constant , 


and the ray remains in the plane dy/dz= constant. How- 
ever, the coordinate system can then be chosen so that 
the ray lies in the z-y plane. Thus, 


dz/ds = 0, n dy/ds = constant , 


d 
o=3 (3), of = cos 0, oY = sin 8, 


so that there follows 
nsin ? = constant. 


Thus, Snell’s law is again obtained. 

From the equations which have been developed here 
—second-order differential equations with given initial con- 
ditions—it is clear that, within the framework of this the- 
ory, the splitting of a ray into two can never occur. If 
a point and an initial direction are specified, then the sub- 
sequent course of the ray is completely determined. How- 
ever, this does not agree with the fact that at the boundary 
between two media, a splitting of the ray always occurs. 
From this, it is seen that Fermat’s principle can only be 
a certain approximation to reality. 


2. THE PRINCIPLES OF MALUS AND HUYGENS. LAWS OF IMAGE 
FORMATION 
We now wish to investigate the properties of the func- 


tion S, where 
Py 


S=|nds=S(P,, Py) . 


Py; 
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From Eq. [1.7] there follows 


d : 
grad 8-82——n (523s), ier. 


d 
grad -32 = (S=-Bs), in 2a 
so that 
grads =F n<* [2.1] 
and. 
|grad S|*= n?. [2.2] 


This last equation is often designated as the eikonal equa- 
tion, and S as the eikonal. These designations are not 
always unique in the literature. 

Instead of the function S, a related function S can be 
introduced in the following manner: Let F be an arbitrar- 
ily specified surface in space, 


= ,(u,v) or F(x)=0 
and let P be a point external to the surface. Then 


-5=5F,P)=|nas, 


where the path of integration between F' and P is to be 
taken along that light ray which is perpendicular to F 
and which passes through P. Our general variational for- 
mulae, Eqs. [1.6] and [1.7], then imply that 


oe ae) Ree ae 


Here, the first term vanishes because the integration path 
is taken along an extremum, while the second term van- 
ishes at the lower limit as a consequence of the boundary 
condition. Thus, for a variation of the point P, there 
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results 


gradpS =n re [233] 


This is Malus’s law. Verbally, it states that if rays have at 
one time been normal to a surface F, then there always 
exist surfaces, namely, those for which S = constant, which 
are normal to the rays. 

There exists an important geometrical relationship be- 
tween S and S which is known as Huygens’s principle: 


The surface S=constant is the envelope of the surfaces 
S(P,, P)=constant associated with each point P, on F. 


Proof. 

Let u, v be the surface parameters. The equation for 
the envelope of the surfaces S= constant is obtained in the 
standard way by differentiating with respect to the para- 
meters: 

oS/du = 0, oS/ov=0, 
so that 


Ox Ox 
gradp,S:~7 = 0, gradp,S-~- = 0, 


from which it follows that just those rays which are per- 
pendicular to the given surfaces strike the envelope. Q.E.D. 

There are two possibilities for describing phenomena in 
an isotropic medium. We have the equations 


and 
|grad S|?=n?, 

whose connection with one another we wish to investigate 
further. It has already been shown that from solutions of 
the ordinary differential equation [1.6], solutions of the 
partial differential equation [2.2] can be constructed. How- 
ever, the inverse of this is also true: if one has a suffi- 
ciently general solution of [2.2], then solutions of [1.6] can 
be obtained from it. 
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We consider a two-parameter family of surfaces 
F(x°, a,,a,) = 0 and the corresponding functions 


SPP) = se, a, Ay) , 


which, because of Eq. [2.3], also satisfy Eq. [2.2]. Upon 
differentiation with respect to a, and a,, there follows 
aus. = ~ Os ; 

a me Va 
that is, the surfaces oS/da,,—constant are perpendicular 
to the surface S=constant. If it is assumed that the two- 
parameter family of surfaces F is so chosen that the 
surfaces 0S/da,=constant and oS/da,—constant do not 
coincide, then they intersect to form a curve which is per- 
pendicular to S=constant and parallel to VS. This curve 
is a solution of the ray equation. This may be seen as fol- 
lows. If we write dS/Oa,=6,, then the equation of such 
a curve becomes 


Dy = fi(8, 1, Gey D1, Bg) . 
Remembering that VS is parallel to dx/ds and using Eq. [2.2], 
we can write VS=ndx/ds, from which follows 


pa 9 oe » dx, 0 0S _ oS a8 
ds a) = “ds 0x, 0", + Oa, 00,02, 


<i I 2) = on 
~ Oa, 2 k 5a, | 7 Oa; 9” "Ba, 


so that 


which is Eq. [1.6]. Q.E.D. 

Our calculation shows us a very important relationship 
between a total differential equation and the partial dif- 
ferential equation associated with it. This connection will 
be considerably generalized later. The relations were first 
discovered by Hamilton. 
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As an application of the variational principles, several 
theorems regarding the exact imaging of surfaces by means 
of extended pencils of rays will be discussed here. 


a. Imaging of a one-dimensional point manifold by means 
of a one-dimensional manifold of directions 


Object-space quantities will be denoted without primes, 
image-space quantities with primes. All rays from P are to 
somehow converge at P’ and, similarly, those from Q are to 
converge at Q'; furthermore, they are to lie in the same plane 
as PO = = dx. If the unit tangent vector to the ray is de- 
noted by #, then, from Eq. [1.7], the variation laws state that 


dS = S(P, P’) — 8(Q, Q') t:dx) —n'(t'-8x’). 
Thus, with |dx|= 51, 
5S = ndl cos (t, 8x) —n’ Sl’ cos (t’, 3x’). 


The right side of this equation must be independent of 
which ray of the pencil is considered. If we denote two 
rays of the pencil by #, and #,, then it follows that 


ndl{cos (t,, 5x) — cos (t,, 5x)}=n'dl' {cos (&,, 8x’) — cos (t,, 8x')}. 
This can be written vectorially as 

n(t, — t,)-3x = n'(t, — t)-8x’, [2.4] 
which is Bruns’s law. In accordance with the derivation, 
the law is, of course, valid only for sufficiently small dx. 
b. Imaging of a two-dimensional point manifold by means of 
a two-dimensional manifold of directions 


The notation is the sameas in a above. In addition, let £ 
be the plane of 5x and t be the projection of ¢ on #. Then, 
tx S< (0% 
and, from a, 
n(t, — T,)° Sx = n' (tT, — T) dx’, 


n(%, — Ts)° 6,4 = n (ey aaa 1) Ox, 
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where, in addition, the t’s have been so chosen that t,—7, 
and t,—t, are not parallel. The same is true for $x and 
5,x. From the familiar identity 


(4x B)-(Cx D) = (A-C)(B-D) — (B-C)(A-D), 


there follows 


(Sx x 5, x) -[(t, — Ti) x (T3 — 7) ] 
=[Sx-(t%— %,)|[8:%° (%3— %1)] — [Six (t.— 7) ][Sa- (ts —7)] - 


Multiplying this by n* and using the identities obtained 
above, we find that the expression 


n?([dxx 3.x]: [(t2— 1) x (%3 — T1)]) 


is equal to the same expression with primed quantities. 
The geometric significance of the first vector product ap- 
pearing in the above expression is that of a surface element 
defined by 5x and 5,x; however, the significance of the 
second vector product is that of the twice-projected surface 
area of the solid angle d2 defined by the rays t,, t,, and t, 
(for sufficiently small At). From this it follows that 


n? df dQ cos (N, t) = n'df' dQ’ cos(N’, t’) , [2.5] 


where WN is the normal to the surface df. This equation is 
known as the Kirchhoff-Clausius law. 


c. Application of the two theorems to rotationally symmetric 
systems 


1. Let surfaces perpendicular to the symmetry axis be 
imaged in the form of similar surfaces by means of a pencil 
of rays which has a finite opening angle 2u and which 
contains the axial ray. 

Consider a section through the axis. Writing y instead 
of 3x and applying Eq. [2.4] with ¢, and t, as the extreme 
rays of the pencil, we obtain the important sine condition: 


ny sinu=n'y’ sinu’. [2.6] 
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This condition is always satisfied whenever the geometrical 
assumption which was made applies. The sine condition 
can also be derived from Kirchhoff’s law of Eq. [2.5]. 

2. Let sections of the axis be imaged by a pencil of open- 
ing angle 2u. If we take t, as the outermost ray of the 
pencil and #, as the axial ray, then it follows from Kq. [2.4] 
that 


nl(1 — cos u) = n'l'(1 — cos wu’). [2.7] 


In general, conditions [2.6] and [2.7] cannot both be sat- 
isfied. This would be possible only for u=u', ny=n'y' 
and nl=vn'l’. In practice (for lenses) this is not the case. 
However, both conditions can be approximately satisfied 
when uw and w’ are small. 


d. Imaging of a surface by means of rays within it, where the 
lines connecting exactly imaged points are the rays themselves 


Let P,, P,, and P, be exactly imaged points. From 


Figure 2.1 


Eq. [2.7], it follows that 
nl,.(1 — cos u) = n'li,(1 — cos uw’) . 
On the other hand, we also have 


nl,,(1 — cos wu) = n'l,(1 — cos u’) 
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and 
nl.3(1 — cos v) = n/U,(1 — cos v’) , 
so that 


Lie lag ites = Lis ths cls ‘ 
Thus, the triangles are similar and their corresponding 
angles are, therefore, equal: 


Wane v9’, nal - 


3. HAMILTON’S THEORY 


a. Special homogeneous form (isotropic bodies) 


In geometrical optics, we have become acquainted with 
the fundamental extremal principle of Fermat 


P, 
[nds —0. 


Py 


On the other hand, we know that in the mechanics of a 
mass point the trajectory is determined by the extremal 
principle 


b) i V2m(E—V(a))ds=0. [3.1] 


Consequently, if we set n?=2m(H—V), then we obtain a 
transformation scheme between the equations of mechanics 
and those of optics such that with each optical problem 
there is associated a mechanical problem. However, the 
analogy does not refer to the time development of the 
phenomena, but, rather, it connects the trajectories of a 
mass point with given initial energy in a potential field 
and those of a light ray in a medium with a given index 
of refraction (as a function of position). 

Thus, in optics as well as in mechanics, the extremal 
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principle of Eq. [1.5] can be expressed in canonical form. 
If we set 


PS ae [3.2] 
then 


With v as a parameter, there then follows 


= s5 = [nas—8[p — de 


=|{p-S+p- <3} ae=|(ap-5 —_—— = E dx) d+: aa 


Now, however, it must be noted that Sx and 5p can no 
longer vary independently. Because of Eq. [3.2], we must 
have 


Py 


sG=0, where G=}(p?—n?*). [3.3] 


Then, using the method of Lagrangian multipliers, it fol- 
lows from 3S=0 that 

dx oG dp oG 

= A Te yh 

dv op dv Ox 


or, after substitution, 


d 


oe dp 
= = Jn — . A 
ap AP? An [3.4] 


Obviously, these are the same equations as [3.2]. By com- 
parison there follows 
2 
ae =j)?n?, from which An=—. 
dv v 


Substituted into Eq. [3.4], this yields the canonical equa- 


tions 
dx oP on 
ead 3.5 
a ae Ox’ | 
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If p is eliminated between these two equations, then the 
familiar ray equation [1.6] results. On the basis of our 
analogy, everything proceeds exactly as in mechanics. If 
we set n? = 2m(H—V), then p becomes the mechanical mo- 
mentum and 8G=0 becomes the law of conservation of 
energy, as can easily be verified by substitution. 

Furthermore, we note that mathematically the three- 
dimensionality is of no consequence in the calculations. In 
fact, we know that the whole formalism can be carried out 
in exactly the same way in the mechanics of NV mass points 
(in which there are 3N dimensions). In optics, however, 
only three dimensions occur. 


b. Most general homogeneous form (erystal optics) 


Furthermore, it is also seen that the specific form of the 
function in 5{nds = 0 is unessential. Thus, in complete 
generality we will set 

P, 
f 
ss = eon [3.6] 
Py 
where uw is an arbitrary parameter, and J is a function 
which satisfies the following conditions: 


~ 


d d 
L=L (a, seey WUptry a seey “te ’ 


du 
i ar ae 
vo Op41 
and 
L(Qiy +++5 Ur Aq, seey Aqs+1) =1L(q:, q:) : [3.7] 


The requirement of Eq. [3.7] is necessary because the 
parameter u is completely arbitrary, and it must be de- 
manded that the variational principle of Eq. [3.6] be in- 
variant with respect to changing the parameter. In the 
case of optics we have f= 2. Euler’s homogeneity relation 
corresponding to Eq. [3.7] is 


, ob 
haga L. [3.8] 
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Again, in analogy to Eq. [3.2], we set 


OL 
= [3.9] 


From Eqs. [3.8] and [3.9] there then follows 


L = Dee: , C35 — 5| (S peat du. [3.10] 


Because LZ is homogeneous of the first degree in the q, 
(according to Eq. [3.7]), then the p, are homogeneous in 
the zeroth degree. If the first f equations of [3.9] are solved 
for qg, (k=1, 2, 3,..., f) and these then substituted into the 
(f+1)th equation of [3.9], then there follows a relation of 
the form 


Prta == DrtilQiy --+y Uptry Dry oory Pry 94-41) ‘ 


Because of the homogeneity (of zeroth order) of p,,, in q, baie 
it follows that p},, cannot depend upon q),,. Thus, anal- 
ogous to Eq. [3.3], there must exist a relationship 


G (Pry oo) Detay hy o2%) Uti) = 0. [3.11] 


The variational principle [3.10], together with the auxil- 
iary condition [3.11], determines the course of the ray. 
The variation becomes 


d 
oe Pie yo du 
dq: 
= Ss op, — aaa eee P* Bas du + Pr ele [3.12] 
k du Py 
With the auxiliary condition [3.11] there finally follows 
o@ F oG 
= [iia =—=— A eet | es 
Oa ae a [3.13] 


On the other hand, it follows quite generally from Eq. [3.6] 
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that 

d oL ob 

du dq. oq 
(Here, there is no auxiliary condition because we have 
worked with the original number of variables without in- 
troducing the p,’s.) Thus, the last equation of [3.13] can 
also be written as 


ee [3.14] 


Equations [3.13] and [3.14] are called the canonical equa- 
tions (analogous to Eqs. [3.4] and [3.5]). 

Here too we can consider the surfaces S=constant (wave 
surfaces), where 


Ps 


S=(Ldu. 


P; 


If, in the variational principle of Eq. [3.6], P, is held fixed 
while P, is allowed to vary on the surface S= constant, then 
we must have 8S=0 by construction (the integral is taken 
along an extremum). However, from the general variational 
formula [3.12], it then follows that 


OL 
an 84=0,  8¢ on S=constant. [3.15] 


This is the so-called transversality condition. In general, this 
will not be the orthogonality relation [1.7] with 5S=0, 
P, fixed. Thus, the rays are not perpendicular but only 
transverse to the surfaces S= constant (see p. 20). 

A differential equation can easily be constructed for the 
function S. From the general variational formula [3.12], 
it follows that 

os 


Ode = Pr [3.16] 
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and this, substituted into Eq. [3.11], yields 


G ay a68. =) Gill 
aq? Odni 9 V1 ery 11] UL. [3.17] 


If it is now assumed that we have a solution of Eq. [3.17] 
which depends upon f constants a, (o=1,2,3,..., f), then 
the theorem that the equations 


[3.18] 


determine light rays is valid. This is the generalization of 
the relation between a total differential equation and its 
associated partial differential equation which was men- 
tioned earlier. 


Proof. 
From Eq. [3.18], 


dy = Pr(U, Ayy w+05 Az, bh, cong b;) [3.19] 


(where we have tacitly assumed that the system of equa- 
tions [3.18] is solvable for the q,; that is, the matrix 
6*8/Oa,0q. is of rank f). If this is substituted into Eq. [3.16], 
there follows 


os 


aa — P);, — Wx( Ay, weey Bry Dy, ones by, U). [3.20] 
": 


Equation [3.17], considered as a function of the q’s and a’s, 
when differentiated with respect to a, yields 


oG oS 


eG 3.21 
Dam ond, fe" 


On the other hand, it follows from Eq. [3.18] differentiated 
with respect to wu that 


O78 dqe _ 


3.22 
— Ba, 0q, du [5.22] 
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These can be simultaneously satisfied only if 


das _ 24 


= : eles 
du Ope [ “J 


Similarly, it follows from Eq. [3.17] when differentiated 
with respect to q that 


oG OG os 


Ode piixead Opi 09109: 
Substituting into dp/du (from Eq. [3.20]), we obtain 
dp, oS dq, oS oG oG 
uu 2 a) ; 

t 


eG, Au 2 &q.0q, Op, Bae 
However, the solutions of Eq. [3.17] given by Eqs. [3.18], 
[3.19], and [3.20] satisfy the canonical equations [3.13], 
which, as shown previously, are equivalent to the ray 
equation. Q.E.D. 
Instead of S, we can here again consider a related func- 
tion defined as follows: 


[3.13] 


i 
S(P) =| Ldu, 
Po 


where P is given and P, is to lie on a hypersurface F(q)= 
constant. The path of integration is to be taken along that 
ray which passes through P and which intersects the surface 
F=constant transversally (condition [3.15]). Huygens’s 
principle again follows, in which S= constant is the envelope 
of the surfaces S= constant corresponding to the points of 
the given surface. (The proof is analogous to that used for 
Huygens’s principle in Section 2.) 

The theory given in this subsection is to be applied to 
light propagation in crystals (see Section 13). It should be 
noted that in this case, ‘‘wave surfaces’? (S=constant) can 
still be defined. Although the ‘‘wave normal”’ p = {p,, ...} is 
indeed normal to this surface, it is in general not parallel 
to the ray direction q'= {q, ...}. 
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c. Inhomogeneous form of Hamilton’s theory 
In the variational principle of Eq. [3.6], the coordinate 


%,+, 18 introduced as the parameter uw. Then, 


dq. a = dq. Wa 
du ’ Gis a 


[3.23] 


The indices have the following ranges: «=1, 2, 8, ..., f; 
i,k=1, 2,3, ...,f+1. Since L(q’, q) is homogeneous of the 
first degree in qg’, there follows 


L= La, veey Vyt1y di. sesy 1) Up-+1 ’ {3.24] 
Ldu =— Ldgs., . [3.25] 
Then, 
oL OL 
Po Bat = Bie’ 


and there follows (from [3.23] and [3.24]) 


oL = Peele, 
Saget ag 
H= > aP.—L, Pri tH = 0. [3.26] 


The canonical equations [3.13] can be rewritten by intro- 
ducing the function H. The second equation of [3.26] plays 
the role of the auxiliary condition [3.11], so that it be- 
comes possible to eliminate the parameter A from Eq. [3.13]: 


G = Pyti t+ A (pay Qa} U1) = 9. [3.27] 
From Eq. [3.13] 
dds 4 oH dP. 408 
du -—S ss Op,’ du Oda” 
ddy+1 ay) AD s+1 ee 3 oH 
du , du Ody+1 


Eliminating A, we get the inhomogeneous form of the 


22 GEOMETRICAL OPTICS Chap. 1 


canonical equations: 
oH ; oH 


da = ae 32 
Va 3 Da ’ Pa Odn [ 8] 
Furthermore, it follows that - 
oH dH 


aaa 3.29 
Odrt1 Ay +1 [ ] 


The connection between the total differential equation, 
the ray equation 
eee 
drt si) Oa." 
and its associated partial differential equation, is also ob- 
tained in the inhomogeneous case, and indeed by again 
specializing the auxiliary condition [3.11] to that of [3.27]. 
Substituting p, = 0S/oq, into Eq. [3.27], there follows 


os os 
sn t Bla Qa) ays) =0. 
Thus, S=S(q,,a,) can be determined, and 
os 
Sa, 


are again the equations for the rays. Of course, upon sub- 
stitution in this manner, H will also be dependent upon 
the f constants a,. 


Isotropic optical systems. The inhomogeneous form is prac- 
tical for application to these systems (f = 2): 


DL=n(i+ a+ #2), or Lda, = nds, 
OL ; moe 
Pa = Te = Mg(L + @1 + #3)t = nda,/ds; 


hence 
n’ — pi — p; = n*/(1+ &i + 23) 
(which can be verified by substitution), and 
L= n(n? — pi — pi, 


i= —- Vi = — 2 


3° 
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We introduce the characteristic function 8, 


os az, 
pee a ag (kK=1, 2, 3), 


and perform the transformation 
V(p*, p) = S— p (x—a) + p*(x*—a’), 
where a and a are arbitrary [A-1]! and 
dS = p-dx— p°®-dx°. 
Then, 
dV =p-dx — p®-dx®— p-dx+ p®-dx° — (x—a)-dp 
+ (x° — a°)-dp® = — (x — a)-dp + (x° — a°)-dp® 


and 
OV OV 
inns ae ene 


We go back to the inhomogeneous form 
P= Vn — pi—-ps=——H. 
Writing 


W(pi, Dos Pry Pa) = 
=V (pi, p2, Vn pe— po, pi, Pe, Vn®— pi — pi) , 


we have 
dW =—>X.dpa+ > Xrdpe, 
with 
Ops Da 
Xp =e — OP ome (%3— G3)= %a— Ag — eat ae (%3— 43), 


and the analogous relation for X3. 


Azially symmetric systems (rotational symmetry). Let the 
z axis be the axis of rotation. We set 


i=, dz/ds = u 
Le =Y, dy/ds = v (direction cosines) . 
ae, dz/ds = w 


Comments [A-1]-[A-7] appear in the Appendix on pp. 151-152. 
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The arbitrary points a° and a are chosen along the z axis. 
According to our general formulas for isotropic bodies, 


ow ow 
Sue ) apes 
with 
A= =— CO ing 
— aV1— a — 08 7 dz 


This shows that (X, Y) and (X°, Y°) represent the co- 
ordinates of the points of intersection of the ray with 
planes perpendicular to the axis at a and a’, respectively. 

The geometric significance of the functions W and V fol- 
lows directly from the analytic definition 


W = S— nt-(x— a) + °F: (x°— a?) 


(where t is a unit vector parallel to the ray). In Fig. 3.1, 


Figure 3.1 
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this is clearly the optical path from N° to N. The func- 
tions V(p°, p) and W/(p%, p$, p,, p,) are called angular 
etkonals. 


4. PHOTOMETRY 


We start from the assumption that for every optical 
process (with the exception of emission and absorption) 
the total amount of light Z is conserved. This quantity 
is also conserved when areas which are bounded by light 
rays are considered. From these premises, Lambert and 
Kepler arrived at the following laws. 

Let F be the radiating surface and f the irradiated sur- 
face. Then dZ= KdQ, where K is a constant dependent 


Figure 4.1 


upon F. If df is an element of the irradiated surface, then 


dO = gpcake. 
v2 
E= —, aL = Edf. [4.1] 


Thus, two surface elements appear equally bright if EF is 


26 GEOMETRICAL OPTICS | Chap. 1 


the same for both; F is called the intensity of illumination 
(or illuminance). To be sure, difficulties are encountered 
in practice when attempting to establish equal brightness 
for different colors.* 

Lambert proposed an additional, although problematical, 
law: the dependence of K on the direction of the light 
emission is 

K=IdF cos, [4.2] 


where I is called the specific light intensity (or photometric 
brightness). Indeed, it can be said that Eq. [4.2] is purely 
formal, since T=I(#). Only for the case of radiation at 
thermal equilibrium is dZ/d? = 0 (isotropic radiator). This 
case corresponds to the fact that a luminous sphere ap- 
pears to the eye as a luminous dise of uniform brightness. 
The combination of both of these laws (Eqs. [4.1] and 

[4.2]) yields 
dL = IdfdF cos « cos #/r? , [4.3] 


dQ = df cos a/r? = sin Pdbdy. 
When integrated for constant I, Eq. [4.3] yields the fol- 


lowing for the total Z radiated into a cone whose opening 
half-angle is U: 


L = |LdafdF cos « cos B/r? 
71/2 276 


=14F| Jeostsind aaay =a Obisin: CU « 
0.60 


For U a2, 
Laat . 


Exact optical mapping. From the preceding discussion, 
the amount of light passing through corresponding surface 


* See J. H. C. MULLER and C. S. M. PovuILiet, Lehrbuch der Physik (Vieweg 
Verlag, Braunschweig, 1925-1934), 11th ed., vols. 1-5. 
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elements must be the same: L=WJL’. Thus, 
Idf cos «dw = I'df’ cos «'dw’. 
However, from Eq. [2.5] it follows that 


n? df cos «dw = n'*df' cos «' da’, 
so that 
ite 1 in [4.4] 


Chapter 2. 
Theory of Interference and Diffraction 


5. GENERAL KINEMATICS OF WAVES 


There exist optical phenomena which contradict the geo- 
metrical ray theory presented so far, and which, in addition, 
cannot be explained by means of a corpuscular theory of 
light. Such phenomena always appear where the index of 
refraction varies very rapidly (but perhaps continuously). 
Among these is the previously considered case of the split- 
ting of a ray into refracted and reflected rays, a phenom- 
enon which cannot be explained, however, by mathematic- 
al epsilontics of surface discontinuities. In addition, the 
diffraction and interference etiects familiar from experiment- 
al physics belong to this category of phenomena. 

All of these phenomena can be treated with the help of 
a wave concept in which the transversality that is specific 
to light will be ignored for the time being. The quanti- 
tative formulation of this theory originated with Young, 
Fresnel, and Euler. 


a. Oscillations 
The simplest form of a harmonic oscillation is 
x= Acos(wit+ 6), 


where w= 2zy is the circular frequency, and v is the num- 
ber of oscillations per second. The principle of superposi- 


28 
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tion holds: 
A, cos (wt + 6,) + A, cos (wt + 6,) = A cos (wt + 0), 
with the relations 
A, cosd, + A, cosé, = A cosd, 
A, sind, + A, sind, = A sind, 


from which A and 6 can be determined. The formulas are 
simpler in complex notation. We have 


a=Ae’ (the complex amplitude) , 
xv = A cos(wt + 6) = Re Aettto — Re aeit, [5.1] 
and the superposition principle, 
a=a4,+4,, [5.2] 


for the superposition of two oscillations. A geometrical 
illustration of superposition is given in Fig. 5.1. As long 


Figure 5.1 


as one deals only with linear operations, everything pro- 
ceeds very simply when written in complex form. However, 
care must be exercised when using quadratic expressions; 
thus, in general, the intensity J is 


I~ (Re ae")? — (Re u)? == xe . with Vhs ae” ; 
Then, 
I~ f(ut¥+u)? = wu + 3 (w+ ut) =| al? + Re ater, 


I~|a|?{1 + cos (2wt + 26)} . 
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Within the bracket the first term is constant, while the 
second term is oscillatory. This oscillating term in the 
intensity cannot be simply ignored, except for the case 
when w is large and a time-averaged value of I is con- 
sidered. 


Vectorial oscillations. Here, uw is a vector instead of a 
scalar: 
x= Reaet = Reu, 


and a = {a,, a2, a3} is a complex vector, i.e., the a; are com- 
plex. The trajectory of a point performing such an oscil- 
lation lies in a fixed plane and is an ellipse (which may 
have the degenerate forms of a circle or a straight line). 


Proof. a. The trajectory lies in a plane. 


dx 
= -— = Reiwae 
dt A 


xxv= Re 3 (ac*t+ ate!) x imac’! = Re= 3 ia*xa, 


since axa=0. However, ia*xa=c is real, since c*=e. 
Thus xxv = we/2 is constant, so that as a consequence the 
point lies in a-plane perpendicular to e throughout the 
entire motion. 


The trajectory is an ellipse. The plane of the motion is 
chosen as the X-Y plane. Then, 
X = 4 (a, et - a, ett) , 
Y= 3 (a, eiwt + as a 


With 
G@=a,+ia,, a*=at—ia;z, 


b=a,—ia,, b* =a*+ ial, 

it follows that 
XA4Y = 4 (ae + b* e-iot) , 
X —4Y = } (det ak e-iot) , 
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We now rotate the coordinate system (X, Y) into (&, ): 


E+ in =(X+1Y)e* = f(a’ ett h’* e-tat) , 
E — ty = (X —1Y)e- = 4 (b' et + g/* et) , 
where 
@ = ae“, a> — ae“, 
bebe , b'* = b*e*, 
and determine the angle of rotation « in such a manner 
that the quotient a’/b’ becomes real and positive. (The 
coordinates corresponding to this rotation are called nor- 
mal coordinates.) 

The special case in which b=0 yields X¥+i¥Y=iace™, 
which is clearly circular motion, and the same is true for 
a=0 [A-2]. Excluding this special case, it follows that 

bja—e--ba—op (real and positive) , 
b/a = oe***, 
whereupon 2« is determined to within modulo 2z, and a, 
therefore, to within modulo z. There follows 


E+ in = 4(a'el*+ pa’*e-), 
é — in — 4 (oa' e+ @’* e~ tor) : 
and, with a’=2Ce, 
E = O(1+ 0) cos (wt + 4) 
n = O(1—¢) sin (wt + 6) J’ 
which are the equations of an ellipse, since eliminating the 
time factors yields 
& i 


[cater (od—or 


An oscillation along a straight line follows for ge=1. Q.E.D. 
Fourier’s theorem permits the representation of periodic 
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processes as a superposition of harmonic processes: 
ft+T) = fo 


+a 
f(t) = > a, Etrot 
tot? 
@ 
,=— t) e inat dt 
a= 2 [fe 


to 


’ [5.3] 


where w=2z/T is the circular frequency. For the for- 
mulas of Eq. [5.3] to be valid, f need not be continuous 
but only to be of bounded variation. 

If f is real, then there follows 


.=O%,. 
Furthermore, because of the orthogonality of the factors 
e'"' we have 


ty 


+7 
pl) —Sa.emlsar == fis [ear—2 3 asa + Sanat 


1 N 
=7{lfolrar— Slash 


for the coefficients a,. From this follows Bessel’s inequality 


1 +N 
Ff lolat> Slat [5.4] 
—N a 
which for N—co becomes the completeness relation 
1 +e 
lf lrat= Slat. [5.5] 


The Fourier formulas can also be generalized to include 
nonperiodic motions by formally carrying out the following 
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limiting process: 


wo—>0, mn—+>oco, nov (fixed), Toco, 


o=Avy, w>>{dy, a,=wA(r). 


Then, 


it oe 
AO ie | f(t) e-** dt 

co . [5.6] 
f(t) =| A(v) et dy 


The completeness relation becomes Parseval’s formula, 


+a 1 +0 
fiaorear == [ito reae. [5.7] 


For this limiting process to be permissible, all integrals 
must of course exist, and the functions must vanish ‘‘suf- 
ficiently fast’’ at infinity. 


b. Waves 


A plane wave progressing in the 2 direction is repre- 


sented by 
u = Acos(kx —wt+ 0), 


for which 
u(x, t) = u(x, t+ T) [5.8] 
where 
Q7 1 
oO = 277 = = y= E 
[5.9] 
27 
p= Q70ph = ein 9 (= oy 
For 


x = vi-+ constant = < t+ constant , 
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the phase of the cosine, 
g=kx—oat+o, 
remains constant. Thus, 
v=o [5.10] 
is the phase velocity of the wave. 
In general, there exists a relation between w and k: 
wo = f(k), dispersion law . [5.11] 


The special case where w/k=c holds for light in vacuum. 
In the three-dimensional case, we use the vectors x and k: 
u= Acos(k-x—wt+6). [5.8] 


Thus, & is perpendicular to the surfaces of constant phase. 
The wave normal n is related to k by 


kan, with n?*=1. [5.9"] 
The dispersion law becomes 
wo = f(ky, ke, ks). = f(k) . [5.11] 


If the medium is isotropic, then we have simply 
w= f(|kl). 


The principle of linear superposition is a general property 
of the waves considered here. It states that if u, and u, 
are two possible wave states of a medium, then so ts ¢,%, + CU, . 
In order that u, and u, be two possible states, the dis- 
persion law must be satisfied. 

The principle of linear superposition is of very far-reach- 
ing significance in wave theory, and it is also well satisfied 
in optics [A-3]. For mechanical waves, however, it is valid 
only for sufficiently small amplitudes. (See, for example, 
a work of Riemann on explosion waves in air.) Because of 
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the superposition principle, Fourier’s theorem permits an 
arbitrary wave state in a medium to be represented as 
follows: 


u=[f fac, k,, ks) cos[k- x — wt + 6(k)] d2k 


= we face ettk-x—ot) 3h 


where a= Ae”. In other words, if any dispersion law of 

the form [5.11'] is given, and if (for optics) it is demanded 

that the principle of linear superposition be valid, then the 

general solution of the problem is given by formula [5.12]. 
With formula [5.12], the initial-value problem 


, [5-12] 


U(x, 0) = U(2) , 


(Ou ; 
(z) U(X) ’ 


is solvable; this is because 


i Re{ {face ek <3; , 
t= Re{ | (— t)wa(k)e* = d8k = tf | fockyakye*=are, 


from which the complex quantity a is determined by 
means of the Fourier inversion formula: 


5 alk) + @t(— m= (5. a) ff tev eth x daer , 
_~ i) Loot) a) —w (—k)a*(—k)] = & JJ f fit eth = d8y 


Group velocity. We consider a wave packet, that is, a 
superposition of many waves of the form given by Eq. [5.12]: 


u(x, t) = Re if | | a(k) eso 43 = Rew. 


36 THEORY OF INTERFERENCE AND DIFFRACTION | Chap. 2 


Here, it is assumed that a(k) is significantly different from 
zero only in some small region around k. The center of 
the packet is denoted by 

o,wdaV 


lee fueav ' 


[5.13] 


We will see that the energy is proportional to u?, so that 
in Eq. [5.13] we have the ‘“‘center’’ of the energy. We have 


u(x, t)=Tw+w*) with w= | | f alk) eh-=-0b 43}, 


e—tot a(k) = gs || fwomrare, [5.14] 


uz = } (w?+ w**) + dwwr*. 


and 


Thus, uw? contains two types of terms, constants and oscil- 
lating terms of twice the frequency. The oscillating terms 
are neglected (time averaging), after which it follows that 


f vww* dV 


a &;(t) == fowtd . [5.14’] 


The numerator becomes 


emit re) 
Bota V = i dinat ik-x) qs 
eso av fa vw Jove i ak, (et*'*) d3k 


= | deanw* i (i * +4 < a) eitk x0 3, 


me (aay ( fare Se are ef w*ad*k) 


where we have integrated by parts and used Eq. [5.14]. 
With Parseval’s formula (Eq. [5.7]), the denominator be- 
comes 


i ww* dix = (27)8 | aa* dek, 
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whereupon 
_, _ fa*i(Oa/dk,) dk | | (Gw/dk,)a* adek 
a(t) = Z ——__, ——____-___ 
a* a d2k ij a*ad%k 
Thus, 
Oz; ow 
> = xe) [5.15] 


is the velocity of propagation of the center of the wave 
packet. 
The group velocity is defined by 


v=—. [5.16] 


In optics, U has the direction of the ray in geometrical 
optics. In deriving the group velocity, we arbitrarily took 
u? as a measure of the energy. Exactly the same result is 
obtained if one takes instead any quadratic form g(w)u? 
(for example (0u/ot)?, since 0/dt~a). 


Isotropic case. We consider a wave packet, Eq. [5.12]. 
For the isotropic case we must have 


w= wk), k=|k|. 


If the packet is sufficiently small, then 


a o ov os 
vr vu=- , aaa, = ° A 


ov 
oA” 
In the isotropic case, for any arbitrary individual wave 
of the wave packet, we can always set 
Vutku=o0, [5.17] 
since for isotropy the relation 


V2 = — k? 


holds for the integrands in Eq. [5.12]. On the one hand, 
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Eq. [5.17] is special in that it is limited to a process which 
is sinusoidal in time; on the other hand, however, it is 
general, since no special dispersion law has been assumed. 
It is also an expression of the superposition principle, but 
only in so far as it refers to wave states of the same fre- 
quency. 


Special case. Let the medium be free of dispersion: 
»= U=ce for all waves of the packet. It is then true for 
all waves of the packet that 


(l= = = constant . [5.17] 


From formula [5.12] it follows that for each wave of 
the packet 


aa = — wt = — ck. [5.18] 
c 
Thus, upon substitution into Eq. [5.17], the wave equation, 


Vite . [5.18'] 


follows. 

From Eq. [5.12], which is the general expression of the 
principle of linear superposition for wave kinematics, we 
have derived a differential equation by taking as a basis 
the quite special dispersion law of Eq. [5.17’]. In this spe- 
cial case—and only in this case—Kqs. [5.12] and [5.18] are 
equivalent. The linearity of the superposition principle 
[5.12] is reflected in the linearity of the differential equa- 
tion [5.18’]. 


One-dimensional case. We consider a single position co- 
ordinate # for the special case where Eq. [5.18’] is valid, 


Ou 1 eu ™ 


oe OP ee 
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According to Eq. [5.12], the solution is 


+a 
2 = | a(k) elk dk + complex conjugate 


—@O 


= fave etkiz—ch) qi + a(— k) e-thx+e qh 
° 0 + complex conjugate . 


Since a(—k) is completely independent of a(k), 
2u = f(a — et) + g(x + ct) [5.19] 


follows as the most general solution. Equation [5.19] can 
also be derived without Fourier integrals by the direct 
solution of the wave equation [5.18’] for this case. 

In the same way, for spherical waves in an isotropic, 
dispersion-free medium, it follows that 


ru = f(r—ct)+g9(r+et). [5.19] 


To show this, one need only to substitute v=rw in the 
wave equation [5.18], and Eq. [5.18"] (with « =r) follows 
for v. If there is to be no singularity at r=0, then one 
must require 

ru = f(r — et) — f(—r—et). 


6. REFRACTION, REFLECTION, INTERFERENCE 


We now apply the general wave theory treated in Sec- 
tion 5 to optics. For this to be at all permissible, it must 
be shown that the ray theory treated in Chapter 1 is con- 
tained within the wave theory, or that it is at least a cer- 
tain limiting case of the latter. This will now be shown 
for the fundamental laws of refraction and reflection. 

Let a bounding surface F be given which separates regions 
characterized by k, w and k’, w’. Our considerations will 
apply to a plane wave in an isotropic medium. 
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We require continuity of the phase at the boundary F. 
This requires 


Os [6.1] 


Furthermore, k-x =k'-x' for x = x’ on F. Thus, (A—f )-x 
=0 for x on F. Hence, (k—f') must be parallel to the 
surface normal N. This, however, is the law of refraction: 


(kk) x N= —k x — 0; 


so that 
ksind = k’sin?’. [6.1] 


Comparing this with [1.1], we see that 
nin!’ = k/k'=2']A . [623] 


We have remarked that in the ease of isotropy, the phase 
velocity is r=a/k. Substituting this into Eq. [6.2] and 
then using [6.1] gives 
n 
a [6.3} 
The absence of dispersion in vacuum is an experimental 
fact; consequently, in our theory a universal velovity e 
must be assumed for light in vacuum. Furthermore. if 
nm=1 in vacuum by definition, then it follows that fer an 
arbitrary medium with imdex of refraction n (with respect 
to vacuum) we have 


¢ t 
4 [6.3"] 


= 
This is the relation between refractive index and phase 
velocity which can be veritied by direct experiment. The 
numerical value of ¢ must be determined experimentally; 
it is not predicted by the theory. 

Continuity of the wave field at the bounding surface F 
is also required for the reflected ray. Just as before, it 
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follows that w=w’ and (k—k’)-x—0 for x on F. How- 
ever, from the above we now have n—n’' (the same me- 
dium), so that |k|=|k’|. Thus, (k—k’) is perpendicular 
to F, or (t—t’) is perpendicular to F, where |t})=1. This 
is the law of reflection. 

In order to show more generally that the ray theory is 
contained in, or at least represents a limiting case of, the 
wave theory, let us now assume an inhomogeneous isotropic 
medium. If processes harmonic in time are assumed, then 
Eq. [5.17] is valid for such a medium: 


V2u+ k?u=0, [5.17] 


where k is a function of position. Introducing the refractive 

index in accordance with Eq. [6.2], we obtain 
V2u+kon?u=0, [6.4] 

where k, is the value of k corresponding to vacuum. For u 


we assume 
u= Reees-o , [6.5] 


where S=S(x) is complex. In addition, u,=e™* must be 
a solution of Eq. [6.4]. We have 


Oy os os 
= tk : 
sat = {( i) + sal * 
Upon substitution into Eq. [6.4] and simplification, the 
equation 
1 
— (VS)? + is V?27S +n? = 0 [6.6] 
@ 
follows for S. 
In geometrical optics there is a similar equation, namely, 
— (VS,)? + n?.= 0 


(see Eq. [2.2]). Furthermore, in geometrical optics there 
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are also the relations 

VS, = p = nt, Vi5, = divip . 
In order that Eq. [6.6] be equivalent to the equation from 
geometrical optics, we must require that 


— (div p) <p. [6.7] 
0 


In the one-dimensional case, this yields 


—|—|<p? g <k,n’*, 
ko | da : dx 
A dn 
— : 6.8 
Ee Ga | ps 


Thus, the change in n over the distance of a wavelength 
must be small—in that case wave theory and ray theory 
give the same results. Condition [6.8] is certainly satisfied 
if the wavelength 4 is very small. Thus, we can expand S 
in powers of A, or of 1/k, to obtain directly 


1 1 
S=8,+—86,4+..+48,.+.., 
ko ko 


where S, is a Solution of the equation from geometrical 
optics (the first approximation in our expansion, if all terms 
with powers of 1/k, or higher are neglected). 
Accordingly, it is seen that if A—-0 is assumed, then for 
the case of isotropic media our wave theory describes phe- 
nomena in the same way as geometrical optics does. The 
question now arises whether the wave theory is a better 
approximation to nature than is the ray theory for cases 
where the finiteness of A is of experimental significance. 
Interference phenomena show that we are on the right 
path when using the wave picture for light. We set the 
light intensity introduced in Section 4 (the illuminance, 
for example) proportional to the square of the wave func- 
tion. Then, in our new theory, in contrast to geometrical 
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optics, the intensities are not additive; this is, indeed, ob- 
served to be the case. For example, let two (harmonic) 
waves be superimposed at a given point. If 


u, = Rea,e-*, 


U, = Rea,e—* 
are the individual excitations of the two waves, then 
U,+ u, = Re(a,+ a,)e-* 


is the total excitation at the point under consideration. 
Thus, there follows 
I, |a,|*, T,oc|a,|*, 
Ie |a,+a, |* = (a,+ a,)(a; + a3) ’ 
so that 
f=7,+-1,-+ 2V Fi,cosd, [6.9] 


where 6 represents the difference of the arguments of a, 
and a,. If this phase difference is well defined, then one 
speaks of coherent wave trains. In this case, the nonad- 
ditivity of the intensities is evident. If, however, we aver- 
age over all possible 6, that is, light from separate sources, 
then it again follows that J=/,+1,. 

Coherent light can be produced in Fresnel’s mirror ex- 
periment. In Fig. 6.1, Z is a light source which produces 


Figure 6,1 
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the virtual sources L, and Li, by reflections from mirrors 
S, and S,, respectively. We calculate the intensity at the 
point P: 


r= PL,=PS,L, ri=a’+(sd—p)’, 
i ee 7, = a + a p)’, 
rg — 1 = 2dp, 


d 
1+, 2a, so that none, 


Thus, there follows 


__ 2 dp 
es 


which with Eq. [6.9] gives 
I = 21,(1+ cos 6) = 41, cos? 6/2 . 


From this we infer that intensity maxima are obtained for 
6=22n (n integral), that is, for p=and/d. In between, 
for p — (Aa/2d)(2n-+-1), are points at which the intensity 
is zero. 

An additional experiment for producing interference phe- 
nomena is the combining of rays which have the same 


Interference 


Plate 


nh 


Interference 


Figure 6.2 


Sec. 6 | | REFRACTION, REFLECTION, INTERFERENCE 45 


inclination after having been reflected from parallel sur- 
faces (Fig. 6.2). A calculation yields the following result: 
The condition 

2d cos y¥ = mj! 


implies maxima below (transmission) and minima above 
(reflection), while for 


2d cos y = (m + 4)2’ 


we have minima below and maxima above. 

Newton’s rings are another interference phenomenon. 
This experiment demonstrates that the wavelength de- 
creases as the color of light changes from red to violet. 


7. THEORY OF DIFFRACTION 


a. General solution of Kirchhoff 


According to geometrical optics, small openings should 
always cast sharp shadows. Experimentally, however, the 
familiar phenomena designated as diffraction are observed. 
Thus, we attempt to explain these phenomena on the basis 
of the wave concept developed here. 

So as not to have to calculate with the general formula 
[5.12], we assume that the entire wave field is harmonic 
in time, with the time factor e~“’. It is clear (because of 
continuity) that if at one point of the wave field we assume 
an excitation of the form e“* unbounded in time, then 
this must also be done for all other points, and with the 
same w. In this way, a mathematically correct formulation 
is obtained; in nature, however, this idealization is only 
approximately satisfied, since in reality one must always 
deal with wave packets (of finite duration). 

With this assumption, upon requiring an isotropic me- 
dium, Eq. [5.17] holds: 


Viu+ktu=o0. 
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In diffraction phenomena, the screens have the significance 
of boundary conditions for w and ou/dn. Equation [5.17] 
is to be solved for these boundary conditions. It should 
be mentioned at the outset that this mathematical prob- 
lem is exactly soluble only for very special cases. For the 
most part, we must restrict ourselves to approximations. 
Here, we will treat the approximation of Kirchhoff. 
We have Green’s theorem for a closed surface: 


Ou Ov 
fo vu veyar =p (o5t— us) af. 


Va FQ 


In this formula we substitute a solution of [5.17] for wu, 


and for v we set 


(= 


Tp , 
(The light source Z is outside the surface F.) The point 
P is an arbitrary field point and Q is the point which varies 
over the region of integration. Thus, 

Vy»+kv=0, 
so that for Q+P 


- uV%—v V2u=0. 


If P is within the region of integration, then it must be 
excluded by means of a small sphere: 


) ) 
bos — us] adf=0 for P outside V, 
F 


@ @ 
bose —u) ay = 0 for P inside V. 


F4-sphere 


In calculating the contribution of the small sphere, we note 
that 6/on = —0/or on the sphere. For r—+0, only the term 
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is significant. On the sphere, df =7r?dQ, and there follows 


Ou ov 
fat (ost 0 52) = — Gd Quy = — trp. 


sphere 


Thus 


Ou ov 
b(0 5 ‘ us) uy 
PF 


4nu, for P inside the surface , 
Peay 


0 for P outside the surface. 


We now let a part of the closed surface go to infinity. 
Furthermore, we require that the solution of Eq. [5.17] 
go asymptotically as 


etkr 


U~m~ 
7 
at infinity, and that there be no contributions of the form 
e **, This is a boundary condition at infinity, the radia- 
tion condition. Thus, physically, only outgoing waves are 
allowed. With this, however, it follows for the integrands 
that 


nee 4, rPdQ—>0 for r>oo. 
on on 

Thus, we need only to integrate over the finite part of 
the surface. There follows 


e. [7.2] 


4nu, if P and L are on different sides of F’, 
0 if P and L£ are on the same side of F’. 


Here, LZ denotes the light source. Formula [7.2] is exact 
if one already has a solution of the problem. In actuality, 
however, one usually knows neither u nor ou/on on such 
a surface; also, these quantities cannot be given in advance. 
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Kirchhoff’s method. The surface is divided into screens 
and apertures: F=.S+A. Then, Kirchhoff assumed that 


ou? ; 

in A: we, ou = ae (undisturbed propagation); 
Ou 7 : wu 

on S: w=0, aaa 0 (vanishing light excitation). 


The quantity u° is to be understood as the form of u ap- 
propriate for undisturbed propagation of light: 


tkr 
‘ — Ee 12a 
u(Q) = ug = constant 


’ 
Tre 


where LZ denotes the light source. 


If Kirchhoff’s hypothesis is substituted into Eq. [7.2], 
then there follows 


1 Cun ig OW 
= w(P) — = |(oF2— ag 2) df. [7.3] 


This is Kirchhoff’s solution of the wave equation. As al- 
ready mentioned, this solution is not exact since the hy- 
pothesis is certainly not correct. If uQ and ouj/on are 
calculated from Eq. [7.3], then the result is not at all 
what was put in. Also, the correct boundary conditions 
are not at all the assumed ones, since on the screens one 
can specify only one of the quantities u or du/dn (not both), 
while at the apertures neither can be given. With this 
hypothesis, Kirchhoff believed that he had found the so- 
lution for a completely ‘‘black’’? screen, but such a screen 
cannot even be defined. 

Kirchhoff’s solution is useful when the apertures are large 
as compared with the wavelength, and when the diffrac- 
tion angle is not too large. 
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If the apertures and screens are interchanged, then Ba- 
binet’s theorem follows from Kirchhoff’s solution: 


u,(P) + u(P) = u(P). [7.4] 


(Here, u, is the solution for the screens at the positions S, 
and u, the solution for the apertures at S.) For example, 
for the case of assumed ideal imaging by a lens system of 
a light source L onto an objective screen (on such a screen 
Up=0 everywhere except at a sharp image point L’), the 
light distribution (~w?) is the same whether a given dia- 
phragm is inserted, or whether another is used which has 
screens where the first had apertures and vice versa. 


Improvement of Kirchhoff’s solution. The function v can 
be chosen somewhat differently. We assume that the dif- 
fracting surface EF is planar, and introduce the point P’, 
which is the image of P in EF. Then, we set 


e* "Pe et? Pa 


Tre Tp’ 


(where Q is the variable point over which we integrate). 
The new v has the properties 


ve=0 
ov 0 et*rre 
on On pq 


on £. The new solution wu, is 
Pe ra eikr 
4nu,p = 4nut -[- U9) 2 - 


8 


df. (75) 


This is related to the old solution by 
“= u(P) — u(P’) + w(P’). 
Another possibility is to set 
era eft Pte 


v1= + 3 
Tpe Tp'e 


50 THEORY OF INTERFERENCE AND DIFFRACTION | Chap. 2 


from which it follows that 


etkr 
v= 
r 
ov 
puree () 
on 
on £, 
0 pikr 
Antip = 4nu® — 2 = = af, [7.57] 


and, correspondingly, 

Up = u(P)+ u(P’) — u(P’). 
However, nothing essential is obtained in this way. The 
advantage of these solutions is that only one boundary 
condition need be specified, instead of the two that had 
to be specified before. 


Consequences of Kirchhoff’s solution. We now substitute 
the following functions uw and »,: 


Vg = efe/r 
Ug = Aetn/r, , 
CE) = Aes 
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However, 


0 fe* d ett ., etkr it 
a (| = cos (n, r) a geen (n, r)ik a (1-z) , 


Thus, from Eq. [7.3] there follows the formula 


A fetkrotry) — 
Up = — | ——— tk 
4a oN 

A 


x{e08(n, r,) (1-5) — cos(n, ry) a : ] 


kr, 


df [7.6] 


for the effect of a point light source. 
We now make the following simplification: 


tie LL, kr 1 (A very small compared with r) . 


In addition, we consider only small diffraction angles—this 
is done in consideration of the experimental fact that at 
small distances away from the geometrical shadow bound- 
ary, the light distribution is that given by geometrical 
optics. Now, because of the smallness of A, the factor 
ert) ig a very rapidly oscillating function, while the other 
remaining factor, cos(m,r,)—cos(n, ro), varies slowly (we 
have already taken 1/ikr~0). Thus, we can assume this 
remaining term to be a constant and replace it by — 2 cos6, 
where 6 is an appropriately chosen average angle. Then, 
from Eq. [7.6] there follows 


etkrgtty) 


up = AS cos 


A 


ToT 


etk(r+ry) 


[7.7] 


TT, 


= w +A cosdl 
27 


s 


If we assume that the diaphragm opening is such that 
a plane can be passed through it, then Eq. [7.7] can be 
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expanded about a null ray. Let M be a point in the open- 


ing and (&,7, 0) the cartesian coordinates of the integra- 
tion point Q. (We have taken M as the origin of our co- 


Figure 7.2 


ordinate system, with the «-y plane in the plane of the 
diaphragm opening.) Then there follows 
m=@—é+y:—nP +e, i=0,1, 
Ria at+yt+at. 
Thus, 
T= Ri—2(@E+ ym) +B +7. 
However, we know that 
a 


Vir & 
ioe oo 
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so that 


aE + Yin rs +p (we +yn) 


I M2es TMs 
, R, OR, Oy lia 


This substituted into Eq. [7.7] yields 


tk — 
ie AS cosé “BJ festemagan [7.8] 
where 
LE+Yon , +7? (HF + yn)? 
@(E,j) = — — te. 
(€,7) R, se oR, OR} +. 
ee eee, 


R, 2R, 2Ri 

With formula [7.8], the light distribution for a diffrac- 
tion phenomenon can be determined (in so far as Kirch- 
hoff’s solution is assumed to be correct) by direct calceu- 
lation. The more terms considered in the expansion of the 
function ® of Eq. [7.8], the more difficult it is to carry 
out the calculation. If only the linear terms in @ are con- 
sidered, then physically this implies that both the light 
source and the field point are infinitely far from the dif- 
fracting aperture. The quantities 7,/R,;=«, and y,/R;=f; 
(t= 0,1) become direction cosines, and Eq. [7.8] becomes 


Up = constant: | Ete Elantadtn(BotPvidedy . [7.8'] 


In practice, one would of course not assume LZ and P 
to be very far away; instead, one would assume that par- 
allel rays had been achieved with lenses. This case is 
known as Fraunhofer diffraction. 

If the quadratic terms are also considered, then we have 
Fresnel diffraction. We will treat some cases of both types 
of diffraction. 
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b. Fraunhofer diffraction 


1. Diffraction by a slit. Equation [7.8] is valid for the 
Fraunhofer case. If the slit (of width a) is infinitely long, 
then we have a plane problem: 


Bothi=9, 

. M14 
% = siny = cos(Z—y), 
a, = —sing = eos(7+ 9), 


pe = k(x, +a) = k(sin yp — sing) . 


Figure 7.3 


With these relations, Eq. [7.8] yields for the slit 


+al2 


constant ; 
mane {en tnale =. efusia} F 


Up = constant: | e~#F dé — 


—al2 
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The intensity is then 


|w|* = constant (“nasi y. 
pe 


T 


Figure 7.4 


If the intensity of the direct ray is called J,, then 


|e|* = 1, Grae [7.9] 


Thus, we have minima for 


a 
Cy me Se ty ess 


or when 
a(siny — sing) =n. 


2. Diffraction from a circular aperture. The circle is placed 
in the é-y plane. Then, Eq. [7.8’] becomes 


u = constant: | e~#eota+nBtAdédn . 


circle 


We now assume that the light is incident perpendicular 
to the diaphragm so that 4,—f,—0 and the 2 axis be- 
comes the symmetry axis for the entire phenomenon. Fur- 
thermore, we introduce the diffraction angle #, and since, 
on symmetry grounds, we see that wu is independent of 
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azimuth, we set {,=0 so that « —sin?. There follows 


wu = constant: i Cane ae dy. 
circle 
Here, we must integrate ¢ over the range —E<&<é and 
ym over the range —a<n<a. If we perform the first in- 
tegration, there follows 
7 eee 
n= constant: | sia sin (ké sin #) dy . 

Now, however, £ can be expressed in terms of the radius 
of the circle and a polar angle: 


€é=acosp, 
9 =asinB. 
4 


Figure 7.5 


Upon substitution there follows 


+712 
constant [. : 
= eae (ka sin 8 cos B) cosB dp. 


—n/2 
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We know, however, that the Bessel function of order one 
is defined by 


+2/2 
J, (2) == sin (z cos B) cosBdB . 


—n/2 


Thus, for «w we get 


u = constant Pultwsiny) ; 
sin} 
and for the intensity 
4 ' 
Io |u|? = constant galas) 
sin? # 


For small z, J,(z)/e~1, so that if the intensity of the un- 
diffracted ray is denoted by J,, we have the following: 


___ (dil(ka sin 8)\? 


The first zero of J,(z) is at z=1.227. Consequently, the 
first dark ring occurs at the diffraction angle 3,, where 


sind, = ——- = 0.614, [7.10] 


Asymptotically, the Nth zero of J,(z) is at 2y~2(N+4), so 
that for the minimum it follows that 


sin dy~(W +3) a 


2a 
3. Diffraction from a grating. By a grating we mean a Sys- 
tem of m (infinitely narrow) equally spaced, adjacent slits. 
Thus, the integral in Eq. [7.8’] goes over into the follow- 
ing sum: 


m1 
u= constant > e-tdkartan | 
p=0 
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or, if we again set 

p= k(%+a,) = k(siny—sin¢) 
(see Fig. 7.3), 


m—1 Ener e-impt —_ | 
u= BOL pet — constant ra 
Thus, 
1— cos 
Toc |u|? = constant ones 
1— cospud 

; 9\2 

= constant (See [7.11] 
sin ud/2 


The principal maxima of order N occur when the numer- 
ator and denominator both vanish: 


= Nn, (ata)o=N. [7.11] 


At principal maxima I~m?. The intensity vanishes for 


d 
a Ne and ee 


2 2 


#0. 


From this, we conclude that the shape of the diffraction 
pattern is as shown in Fig. 7.6. If the width of the slits 


. 


Figure 7.6 


in the grating is taken into consideration, then we obtain 
the product of the formulas for slit and grating. Gratings 
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are the most important instruments for spectral analysis. 
They permit the measurement of wavelengths. 

Technically, two neighboring wavelengths are considered 
to be ‘‘resolved”’ if the principal maximum of one falls 
exactly on the first minimum of the other. 

According to Eq. [7.11], (sing—siny)(d/A)=N, so that 
for normal incidence the principal maximum of Nth order 
occurs at the diffraction angle 


: A 
sin yp = in ; 
and the first zero after that is separated from it by the 
angular interval 


: A 
A(siny) = aii 


On the other hand, if the wavelength is varied then the 
center of the principal maximum is displaced by 


A'(sin y) “4 N. 


If we are to have A’=A, then 


ma 


AQ. 


A= 


Thus, the resolving power A =A/AA/ is equal to the prod- 
uct of the order times the number of lines in the grating. 


4. The resolving power of a microscope. The Fraunhofer 
theory shows how the concept of resolving power can be 
defined. We consider the example of the microscope. 

Luminous object. The point P can no longer be differen- 
tiated from Q if P is so close to Q that P’ falls within the 
principal maximum of the diffraction pattern of Q. 
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If the object plane # is exactly imaged at the image 
plane Z’, then the sine condition 


Insind = ln’ sind! ~ In'd’ 


Objective 


|.~____——-b 


Object plane Image of the Image plane 
entrance pupil 


Figure 7.7 


must be satisfied. In order that P be differentiable from Q, 
with the image of the entrance pupil as the diffracting 
aperture, one must have 


l aM A 


= = oul = —= 6’ 
OP 0.6 ZA? but > 6’, 
so that 
A 
US0.61—— 
> 6 n'd'’ 


which, with the above-mentioned sine condition, gives 


U'n'd! Ao 
~ nsind se nsind~” 


Illuminated object. A grating with grating constant d 
is assumed as the ‘“‘standard” object. Many orders of the 
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diffraction pattern produced by the grating must be in- 
cluded for an image to be obtained. Let y be the angle 
of the first principal maximum. Thus, two coherent light 
sources are formed at C and C, in the focal plane, and 
these sources produce a system of lines with d’=Ab/A at 


| 
Grating 


Objective 
Object plane Focal plane Image plane 


Figure 7.8 


the image plane EH’ (see Fresnel’s experiment, Fig. 6.1). 
Therefore, we get an image whose structure is similar to 
that of the object only if at least one diffraction maximum 
of the grating is included. Thus, 


p< :(opening angle)=6, 


A : 
sing = ad wind 
Ay 
ie nsind’ 


where nsin6d is called the numerical aperture. A factor 
of 4 is gained when the system is arranged so that the 
direct ray hits the edge of the lens. If one is satisfied with 
only the indication of a particle without wishing to know 
its shape, then it is only the light intensity that deter- 
mines how small the particle may be. 
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The concept of resolving power can also be applied to 
other optical instruments. Thus, the calculation for a step- 


Emergent 
light 


Incident light 


Figure 7.9 


grating (Michelson’s transmission echelon) yields [A-4] 
nm—1 dn 


where m is the number of layers. (For additional examples, 
see the Supplementary Bibliography.) 


c. Fresnel diffraction 


We now assume the light source Z and the field point P 
to be finite distances from the diffracting aperture. Then, 
in Eq. [7.8] terms up to at least the quadratic ones must 
be considered. At the same time, we locate the origin of 
the coordinate system at the point of intersection of the 
connecting line LP with the plane of the diffracting aper- 
ture. (For fixed L we have a different coordinate system 
for each point P.) In this way we achieve the result 
that in Eq. [7.8] the linear terms in ® drop out. If we 
now set ¥,=Y¥=0, then «=—a,=sind, and there remains 


(ek i 


; : {t12] 
u(P) = — ik cond enias ent ag dy 


27k, Rk, 
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We make the following change of variables: 


il 
Katee cos? 6 = 2? | 


1 1 ’ [7.13] 
ee | ak — arn? 
“ (zt z) sey 
so that 
al 1 
k (z + z) cosdd&édy = adudv, 
Rh, 
with which Eq. [7.12] becomes 


4 ethRot+ RD 


T= Saree 


eitrr/2yut auf evan : [7.14] 


The limits of integration in Eq. [7.14] are easily trans- 
formed. If the original limits were -+-co (the case of a 
slit), then in the transformed integral they are also + 0o. 
We have 


+o (0) 
2 
et(a/2)v? q. — 9 et(a/2)u3 du == = etinl4) — 1 4 i 
i v= 2 va ‘: 


—o 0 


The same integral between finite limits, 


8 


| etl2ut Gay, [7.15] 


s—D 
is called Fresnel’s integral and is tabulated. 

Example. Diffraction from the edge of a semi-infinite 
plane. 

The second integral in Eq. [7.14] clearly has the limits 
-+co, while the first integral is to be taken from a finite 
limit to infinity. For the discussion of the integral [7.15], 
we consider the curve defined as follows: 


X =| cos (Fw) : Y= a0 G u) du. [7.16] 
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Since (dX/ds)?+ (dY/ds)?=1, then s is the are length. If 
we set dX/ds=cost, dY/ds=sint, then t is the angle 
the tangent makes with the X axis. Then, however, the 
radius of curvature is 9=ds/dt, and it follows from 


Eq. [7.16] that 
I 
Ce [7.17] 
The curve which is defined by Eq. [7.16] is the familiar 
Cornu spiral. It can be constructed approximately, for 


example, from Eq. [7.17]. (See Fig. 7.10.) 


Figure 7.10 


If we consider the integral of Eq. [7.15], then we see that 
& 2 
few du| = X?+ V2 


[7.18] 


esis au = (X= xX ay aes 


3! 
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In Fig. 7.11, the coordinate system has been chosen so 
that the projection of ZP on the plane falls along the a axis, 
which is perpendicular to the edge of the plane. We thus 


Sereen 


L 


Figure 7.11 


see that the limits for § are —oo and #, so that from [7.13] 
the limits for u are —oco and x cosdfk(R,+R,)/aRR,}?. 
From Eq. [7.18], we recognize that the intensity ratio 27/I, 
[where J, corresponds to undisturbed light propagation, 
e* Roth) 1(R,+-R,)] for u(P) = u(s) is obtained by traveling 
along the Cornu spiral a distance s from the origin (and 
with the proper sign), and then taking the square of the 
distance from the point thus obtained to the end of the 
spiral at F. The case s<0 corresponds to geometrical 
shadow, while s>0 corresponds to the geometrically illu- 
minated side. In particular, for 


S=-+00, (UU = 1 
c—0, Lis (shadow boundary) ; 
io 


s=-—-OO, 


The result is shown in Fig. 7.12. On the shadow side, the 
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intensity decreases monotonically; on the light side, how- 
ever, the intensity oscillates. 
If instead of a semi-infinite plane we have a slit of width d, 


Shadow 


Figure 7.12 


then everything is calculated as before except that the lim- 
its for the second integration are now s—D and s, where 


s = x cosd{k(Ry+ R,)/7R,R,}4 
as before, and 


D = dcosd{k(R,+ R,)/7R,R}? . 


Chapter 3. Maxwell’s Theory 


8. FOUNDATIONS OF THE THEORY 


It turns out that the scalar wave optics developed in the 
previous chapter does not agree with certain experimental 
facts. That is, certain optical phenomena can be demon- 
strated which can be explained only on the basis of the 
concept of transverse vector waves (for example, birefrac- 
tion and polarization by reflection). A light beam is said 
to be polarized if the direction of the vector which will 
now be associated with the propagation of light can be 
specified. Experimentally, it is found that: 

1. Interference between light beams of like polarization 
is the same as for unpolarized light. 

2. Two light beams with polarizations perpendicular to 
one another do not interfere. (From this it follows that 
light has no longitudinal component.) Furthermore, there 
is no interference if the two originally perpendicularly po- 
larized components of a light beam are brought into coin- 
cidence by repolarization. 

3. If uniformly polarized light is repolarized several times, 
then there is interference only if finally the light is again 
uniformly polarized. 

The (experimentally determined) pure transversality of 
light waves and their (experimentally determined) velocity 
of propagation suggest the conjecture that the hypothesized 
oscillating light vector is to be identified with the magnetic 
or electric field vector, H or E. As shown in Electro- 
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dynamics,! these fields also satisfy a wave equation. For 
a charge-free region, and for processes varying in time, the 
Maxwell equations (for ~=1) hold: 


Ea? yg div E=0 
c ot 
10H [8.1] 
a o divH = 0 
¢ ot 
From the familiar vector identity 
V2A = grad div A — curleurlA , 
there follows 
vim See 
c* of 
’ [8.2] 
A eee 
c? Ot? 
from which the phase velocity is seen to be 
c 
— = 8.3 
v= [8.3] 


The successful results which follow from the identification 
of the “‘light vector” with E justify this hypothesis. It 
turns out that what was formerly denoted as the ‘‘light 
vector’? corresponds to the quantity H. For this choice, 
the ‘‘plane of polarization”’ is to be understood as the plane 
which contains H and the wave normal [A-2]. However, most 
of the effects of light are to be attributed to the effects 
of E(see p. 83) since, in accordance with the general electro- 
magnetic law of force, E exerts a force on stationary as well 
as on moving charges, while H acts only on moving charges. 
Thus, we will always operate with E. In order to preserve 
agreement with the previous results, we must set 


n=4/é 


1W. PAULI, Lectures in Physics: Electrodynamics (M.1.T. Press, Cambridge, 
Mass., 1972). 
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in accordance with Eq. [8.3]. In an electromagnetic wave, 
E and H are perpendicular to one another and to the di- 
rection of the ray. 

We now consider a plane wave which satisfies Eq. [8.2] 
and which propagates in the x direction. Thus, all field 
components depend only on the variable (x—vt), so that 


O/oy = d/oz = 0, 0/ot = — ce*0/Oz . 


From this, it is immediately seen that with this hypoth- 
esis two separate solutions of Eqs. [8.1] and [8.2] are ob- 
tained: 


(1) HL, H,A0, 
the remaining field components being zero 


[8.4] 
(2) #,, H,#0, 
the remaining field components being zero 
In these cases, the Maxwell equations [8.1] give 
. a. [8.5] 
(2) H,=+~</ek, 


These two solutions correspond to the two types of polar- 
ization for plane waves; the solutions can be linearly 
superimposed, and the form thereby obtained is still a so- 
lution of Maxwell’s equations. 
As is well known, the law of energy conservation also 
follows from Maxwell’s equations: 
OW _ 


where S=(c/42)(ExH), W=1/8n(ek?+ uH?). 
Maxwell’s equations [8.1] can also be written in integral 
form. In that case, the formulas in Eqs. [8.1] which con- 
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tain the curl are equivalent to 


pH.ds = j {etal | 
} B,d8 = 5 Hat | | 


Equations [8.7] are used to advantage when « becomes 
discontinuous (as, for example, at a surface), in which case 
Eqs. [8.1] break down. In such a case, shown in Fig. 8.1, 


[3.7] 


alg 2' 
6 is 
(a 
di 2 
Figure 8.1 


e has a finite discontinuity at the surface #’, and the parts 
of the integrals over the lengths 6 (6 small) are of the 
same order of magnitude as 6. However, the term on the 
right in the first of Eqs. [8.7] is also of order of magni- 
tude 6, since H, certainly remains finite. Thus, there follows 


2 2 
[Beas | H, as = O(8), 
1 1’ 


and from the limit 6-0 we conclude that the parallel 
component of E (and H) remains continuous across a sur- 
face at which e changes discontinuously: 


E,, H, are continuous across a surface where & changes 
discontinuously . [8.8] 


Up to this point, we have proceeded as if ¢ were a ma- 
terial constant, although this is not the case. This can be 
taken into consideration if the time dependence is assumed 
to be periodic, ~e“*'. Then, 0/dt~—iw, and there fol- 
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lows from Eqs. [8.1] that 


El@ 


—— & = curl 


. [8.9] 
H = eurl E 


This can next be transformed into the form of the time- 
independent wave equation (in analogy to the equation 
V?u-+k2u—=0 of Chapter 2): 


VE+ < E=0 | 
[8.10] 
ViH+  H=0 | 


Whether the system of equations [8.9] or [8.1] is more 
appropriate for a specific case must be decided on the 
basis of the experimental setup. 


9. NONABSORBING MEDIA (FRESNEL’S FORMULAE) 


The system of equations [8.10] possesses solutions of the 
_ form 
E= E, e6* 2-28 [9.1] 


(and analogously for H), if we set 


The phase velocity is then v=c/Ve, and, as was noted 
previously, n= Ve (isotropic medium). 

We consider the case where « changes discontinuously 
across a plane. Since Eqs. [8.10] possess solutions of the 
form [9.1], we choose the initial conditions for the general 
optical wave field in such a manner that they correspond 
to an incident plane wave in the #-z plane, making an 
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angle « with the 2 axis. (See Fig. 9.1.) From Eq. [8.4] 
we have seen that the system of Maxwell’s equations—in 
the case that plane waves may be assumed as solutions *— 


Figure 9,1 


possesses two completely different solutions. Thus, here 
too it will be necessary to treat these two cases separately: 


Case 1. E perpendicular to the plane of incidence, 


Case 2. E parallel to the plane of incidence. 


First, however, we can establish that the continuity of 
the phase at the bounding surface, as shown in Section 6, 
results in the fact that w in one medium is the same as w 
in the other, and furthermore that «—«, and v/esing, 
= Ve,8in%, or sina=nsing,, where n= ~/é/e,. Thus, the 
unit vectors n, n,, and n, of the incident, reflected, and 


* This is to be expected because of the chosen initial condition. 
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refracted waves, respectively, become 
n= {sina, 0, —cosa}, n, = {sing, 0, cosa}, 
n, = {sina,, 0, —cosa,}. 
Case 1. E perpendicular to the plane of incidence 
Form of solution: All fields are multiplied by e**. Fur- 


thermore, k, = kn, andn,-x—=s,. Then, clearly, 
EH, — Aetht ( 
H, = V/é, A cosae » incident wave; 
H, = /¢,A sin ae* 
HL, = Be 
H, = — V/@,Beoosae*" }, reflected wave ; 
H, = Vé,B singe 
oy — COeikss: 
H, = V,C cosa,e**= , refracted wave. 
HH, = V/é,C sin «6% 
For the relation between |E| and |H|, see Eq. [8.5]. On 
the basis of the principle of [8.8], since k,s = k,s,=k,s, for 
2=0, there follows 
from the continuity of #H,: A+B=C, 


from the continuity of H,: (A—B)cos«=— Cncosa, , 


so that 
B sin (a% — &) 
A sin(a+ a) 
C 2cosesin«, 
A ‘sin (a+ a) 


[9.2] 


Case 2. E parallel to the plane of incidence 
We take H in the —y direction, so that E now plays 
the role of H in Case 1 above. 
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Form of solution: (For the relation between |E| and 
|H|, see Eq. [8.5].) 


H, = — Aeths . 

EH, = a," A cosae™?® » incident wave ; 
EL — 6 * Asniaess 

A, =— Be 

E, = — e;1*Beosae™" |, reflected wave ; 
B, = €,1" Bamige*s 

H, = — Cethes 

E, = ex"*C cosa,e*: +. refracted wave . 


E, = €,2C sin a, €**s" 
The continuity of the parallel components yields 
A+B=0C, 
C 
(A — B) cosa = Hy COB Mas 


from which 


__ tan («%— a) 


B 
A tan («+ a) 
C 
A 


2 sin 2a , Brel 


sin 2a + sin 2a, 


Discussion. From Electrodynamics,* the intensity I of an 
electromagnetic wave is proportional to the Poynting vec- 
tor S: 


c 
S= 7, Endl, IT=|S|. 
It is also known that for a plane wave (which is the case 


7 W. PAULI, Hlectrodynamics. 
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here), £ is normal to H. Thus, from Eq. [8.5], it follows 
that for our case 


Toc|E||H|= |H|?/n = n|E}*. [9.4] 


Thus, as long as total reflection does not occur, so that 
the exponential factor e“"°? has magnitude unity, the 
ratios of the reflected and refracted intensities, respectively, 
to that of the incident intensity for the first case are 


) De at Beige) 050 ae y ares Gen = n|0/A|? [9.5] 
from Eq. [9.2], while for the second case 
J eetipey bee = |B/A 9 ys ee = | C/A |2/n [9.6] 


from Eq. [9.3]. 
Of course, at normal incidence both formulas become 
identical: 


j em Bae — {(n — 1)/(n Se 1S » y enon foe = 4n|(n an i? bs [9.7] 


At glancing incidence, «=2/2, it follows in both cases 
that I,,,,:J,,=1. Furthermore, for the second case, it is 


seen from Eq. [9.3] that there exists an angle «, such that 
I .,,:1,,, is zero. This is the case when 


refl+~ inc 


Oy + Ky = 70/2 [9.8] 


(since tanz/2 = oo); this is the case when the refracted 
ray and the (nonexistent) reflected ray would be perpen- 
dicular to one another. This is Brewster’s law. Thus, only 
one component of unpolarized light incident on a reflect- 
ing surface at an angle «,, as determined by Eq. [9.8], 
is reflected. Consequently, the reflected light is linearly 
polarized. 

Since the component of the vector E which oscillates 
normal to the plane of incidence experiences different ef- 
fects than the parallel component, then oblique linearly 
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polarized light remains linearly polarized after refraction 
and reflection, although in general its plane of polarization 
is rotated. 

Somewhat different relations are obtained for the case 
of total reflection. Since formulas [9.2] and [9.3] are purely 
formal solutions of Maxwell’s equations, then they of course 
remain valid for this case too. On the contrary, formulas 
[9.5] and [9.6] for the intensities are no longer valid, since 
the exponential factors are no longer necessarily of mag- 
nitude one. 

In order that total reflection occur, in accordance with 
Chapter 1 we must assume that n<1 and that « is suf- 
ficiently large so that 


sing = nsin a, 
no longer has a real solution for a. Thus, 


a>a,, Where sina,=n. [9.9] 


The angle «, defined in this way is called the limiting angle 
for total reflection. 

As can easily be shown, for sine, to be real and to be 
greater than unity, «, must necessarily be of the form 


= 1/2 —iB, [9.10] 


where f is real (the sign in front of 8 is negative by con- 
vention only). With this substitution, if we completely 
write out the expressions for the field components (inclu- 
ding the exponential factors), we obtain 


Ey — Be-itwt-kisy — Be ~itot—ki2sino—k,z cos &) 


—H, reflected wave, [9.11] 


H 
i [es C e~tmt—kas.) — ( e—iwttik,zcoshB e*sz sinh B 


refracted wave. [9.12] 
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This represents a so-called surface wave, since there is pe- 
riodicity only in the x direction, while in the z direction 
there is an exponential falloff. Upon computation, the 
Poynting vector is seen to be parallel to the plane of re- 
fraction. For both cases of polarization there follows 
Lien tine= 1, 80 that all the light is reflected. On the con- 
trary, it is noticed that for reflection the field components 
for the two cases have different complex factors; that is, 
linearly polarized light oblique to the plane of incidence 
becomes elliptically polarized after reflection. The complex 
factors are just those values for B/A given by Eqs. [9.2] 
and [9.3] when the substitution [9.10] is used. 

With Fresnel’s formulae, a modification of the geometric 
ray theory has been obtained, such that even for the case 
of total reflection a wave in the second medium must be 
assumed, in accordance with Eq. [9.12]. It can be shown 
experimentally that (as is demanded by the theory) light 
is always transmitted through an interrupting layer of 
finite thickness with n<1, even if total reflection should 
occur according to the geometric theory. 


10. ABSORBING MEDIA (OPTICS OF METALS) 


We start once again from the Maxwell equations, except 
that we assume that the medium possesses an electrical 
conductivity. Then, for processes which are periodic with 
e **' the Maxwell equations become 


curl = (—<io +) B 
j [10.1] 
cur] E = ioH 


(since 0/ot~—iw). We can rewrite the first of these equa- 


tions: 


4 
eur) H = (—£ iw +32) g—— Sn" i, 
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where we have defined 


4 
f=. = [10.2] 
w 
Furthermore, if we substitute 
n’'=n(l+ix) (n, x, real), [10.3] 


then upon comparison with Eq. [10.2] there follows 


n2(1— x?) =e | 


[10.4] 
nx = 2na/w 
Thus, Maxwell’s equations of [10.1] become 
1@ 
curl H = — a elE 
[10.5] 


eurl E = CH 


Equations [10.5] are exactly the same as Eqs. [8.1] written 
for periodic processes, except that «’ appears everywhere 
in place of «. Thus, all the previously developed solutions 
of the Maxwell equations are also valid here if only we 
substitute «’ (complex) for ¢ and n’ for n. 

For example, in Eq. [9.1] we had found a solution of 
Maxwell’s equations which, when carried over to the pres- 
ent case, becomes 


E=E,es-0 = with = x = {a,0,0}, [10.6] 


where (in this case) 


Thus, from Eq. [10.6] we get (as a special case of Eq. [8.4]) 


E, — Bi e7(#no! ode et((new/ehe—wt) ‘ [10.7] 


Sec. 10 | ABSORBING MEDIA - 79 
From Eq.[8.5], the corresponding magnetic field intensity is 


BH, a ie e7 nel ox et((nw!cda—ot) ‘ 
where 
He — ae [10.8] 


Equation [10.7] represents a damped plane wave whose 
damping constant is proportional to x. Here too, n (with- 
out the prime) can be interpreted as the refractive index, 
since c/n is the phase velocity. It should be noted that, 
according to Eq. [10.8], E and H no longer have the same 
phase, since n’ is complex. 

If in the above equations we again set o—0, then it 
follows that x—0 and that n=n’' is real. For the case 
of metals, however, o+0. Indeed, o is so large that n2x> ¢; 
that is, x~1 and n is large. 

We now consider the case of a linearly polarized plane 
wave incident on a plane metal surface, with the angle 
of incidence being a and the angle of refraction «,. (See 
Fig. 9.1.) The Fresnel formulae [9.1] and [9.2] for both 
cases of polarization can, in accordance with the above 
remarks, be taken over if « is replaced by e«’ and n by n’. 
Because of the assumptions for metals, it follows that 
cosa,~ 1. 


Case 1. E perpendicular to the plane of incidence 


C 2cosasina, 2 cosa 


an ~ - [10.9] 
A sin(«+ «,) n 


Thus, for the refracted wave, 


a —— exp fi[ky(# sine — nz)—ot] + nk,x2} 


H, =, 2A cosa exp {i[k,(asinw — nz) —at] 
+ nuke} {- [10.10] 
_ 2A : : ' 
=) 4/ 6, =F cose sina exp {i[k,(a sin a— nz) — wt] 


+ nuk,2} 
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For n'->0oo, only H, remains finite (even for z= 0). 


Case 2. E parallel to the plane of incidence 


ee 2 [10.11] 
A sin2«-+sin 2a, 

Thus, 

H, = — 2A exp {i[k,(w sin a — nz) — wt] + nxk,2} 

E, = €°? = exp {i[k,(w sin «— nz) — wt] + nxk,2} [10.12] 


E, = e;? sin« exp {i[k,(x sin « — nz)— wt] + nxk,z} 


for the refracted wave. Consequently, in the limiting case 
of a perfect metal (n'’— oo) the electric field intensities go 
to zero (for 2=0 as well). However, it is clear that this 
must be so, for the current is proportional to E, and we 
have assumed that the conductivity o is infinite. 
Normal incidence. From Eq. [9.7], there follows 


B 


A 


2 


ial 
n+1 


2 (m — 1)? + n? 32? 


Let —s 


die 


(a 1) + mt” 


From this follows the important formula for metals 


7 Es 4n th 2 0/2 
In (n+1)2?+n?x? Qn? n 220° 


In accordance with Eq. [10.4], for «1 there follows 


Lrent € 
1— == |i, | : 
/ i [10.13] 


inc 


Ellipticity of the reflected light. We assume that linearly 
polarized light, with field amplitudes A, and A,, is inei- 
dent on the metal surface. Then, from Eqs. [9.2] and [9.3] 
there follows 


(3) sin (a — @,) (3) tan (« — a.) 


~~ gin(a+ a)’ A},  tan(«-+ a)’ 
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Thus, 


B, ae ais COS (& — a) 
B, Ay cos(a+a,)° L104) 


Since the quotient of the cosines in Eq. [10.14] is, in gen- 
eral, not real, the reflected light is elliptically polarized. 
We consider the special case where A,/A,=1; that is, the 
plane of incidence and the plane of polarization make an 
angle of 45°. Furthermore, we set 


e — Ger, [10.15] 
We compute 
1+ee4  _—_ sina sina, 
1— oe*4 COSA COS&, 
However, 
SI, — -, sing, COS @ = ca 0,16] 


which, upon substitution, gives 


1+ oe4 sin? 
1—oe4 cosa 


(n'? —sin? ~)-*, [10.17] 


From this we see that for 
a= 0: pt AS x; o=1, 
e== 7/2" oe |, AO, ae 
there is no ellipticity. 
The angle of incidence «,, for which 4—z/2 so that 
e'4—j is called the principal angle of incidence. For this 


case, if we neglect sin?«, compared with n’’, there follows 
from Eq. [10.17] 


LA 45 _ Sith ee 1 [10.18] 


1—i1o,, COS a, 0! 


The quantities g,, and «, are measured experimentally, 
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so that the optical constants follow from Hq. [10.18]: 


sin? @,, (; — see 
’ 


n'=n(l bin) = 


~~ cosa, \1+ io, 
n= sin? Xpr (; = 4 
~~ G08 &, \1-+ 8, 
. z | [10.19] 
paar = tan2y, Op =—tany 


It turns out that the ‘‘constants”’ are frequency depend- 
ent. They all approach their static values if w becomes 
sufficiently small. 


11. STANDING WAVES 


We consider a plane wave incident on a perfect mirror 
(H, vanishes for z=0). We obtain the following for the 
two cases of polarization: 


a. E perpendicular to the plane of incidence: B/A=—1 
The superposition of the incident and reflected light waves 
yields 
E, — Ae ivttikte sln&é—gco3 a) ___ A emimttina sin &+2coac) 
= — 21 Ae tottikesine sin (kz cos «) 
H, = V/é, 2Ae- i+ ike sine eos x COS (kz COS &) j _ 
H, = — Vé, 21 Ae itt ikesine gin g sin (kz COs &) 


b. E parallel to the plane of incidence: B/A=1 


H,, = — 2 Ae-twttikzoine ogg (kz cos «) | 
Ei, = — e7'? 21 Ae iatt tkesing egg gy sin (kz cosa) |. [11.2] 
EB, = e718 2Ae—ttttkesina sin & cos (kz cos a) 
Normal incidence. Both cases are the same: 
BH, = — 2|A|sin(wt + 6) sin ke 
H, = V/é,2|A|cos(wt + 4) cos kz Lhe 
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Thus 

1. E has nodes at ke=aN, 2=NA/2. (The first node is at 
the mirror surface.) 

2. H has nodes at ke=n(N+4), e=4(N+4)A. 

If such a situation is photographed, the first node is found 
to be at the mirror. Thus, in photometry it is |E|? that 
is measured. This can be understood on the basis of the 
theory of electrons. 


Arbitrary incidence. For Case a, from Eq. [11.1] it fol- 
lows that nodes occur for kecosa=aN or at z= 4NA cose. 
For Case b, it follows from Eq. [11.2] that for «40 and 
a-~n/2, |E\?=|E,|?+ |H£,|?40, and thus there are no nodes 
present. From this it can be determined that the plane 
previously called the ‘‘plane of polarization” is parallel to 
the magnetic vector H (see p. 68). The vector E is perpen- 
dicular to this plane. 


Chapter 4. Crystal Optics 


12, RELATIONS FOR THE WAVE NORMAL 


In discussing Maxwell’s theory, we have assumed that 
the medium was isotropic. Now, however, we wish to as- 
sume that E is no longer parallel to D. Thus, we must 
write the Maxwell equations in their most general form, 
where, however, we once again assume a permeability w=1. 
It turns out, of course, that at the high frequencies which 
are applicable to light, u(w) is close to one. Thus, we have 


cul =~, div D=0 
a [12.1] 
curl E = —— H, div Hi =} 
From Eqs. [12.1] there follows 
: Anse 
V°E — grad div E = re ; [12.2] 


where divE no longer vanishes. For the relation between 
D and E we must assume 


D=e:-E or D;=)> én, [12.3] 
k 


where e€ is the dielectric tensor. The Poynting vector 
becomes 


Cc 


84 
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and, from the law of conservation of energy, there follows 
(cf. Eq. [8.6]) 
OW 


is c 
ia —divS= — 4g, (H curlE — E-curlH), 


where W is the energy density of the field. From this, 
with Eqs. [12.1], 


eh algers) oD 
“Ot. 8x ot + 7E Ot” peel 


This must, of course, be an exact differential, so that we 
must postulate 


oD OE 
E: a =D~) [12.5] 
for then 
oD 10 (E-D). 


= “Ot 8x ot 
If Eq. [12.3] is substituted into Eq. [12.5], there results 

Yenk EF, =a Dd em Bm Ei , 

tk lm 


or 


> (e% — &x:) EB; EB, =0, 
th 
Esk = Ex: + [12.6] 


Thus, the law of conservation of energy requires that the 
tensor € be symmetric. However, such a tensor can al- 
ways be transformed to principal axes. The characteristic 
equation for the eigenvalues is 


&u— A E12 €13 
E91 Eo — A Ea3 — 0 * [12.7] 
€31 Es €33 — A 


Furthermore, we know that the directions of the principal 
axes corresponding to two different roots of the charac- 
teristic equation are perpendicular to one another. If two 
roots are the same, then the corresponding directions of 
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the principal axes are undetermined in a plane, so that 
they can be chosen orthogonal. For a (real) symmetric ten- 
sor, all eigenvalues are real. 

In what follows, we assume that the transformation to 
principal axes has been made at each point of the crystal, 
so that the tensor € is diagonal and can be written as 
follows: 


&y 0 0 
(ex) = | 0 bs 0}. [12.8] 
0 0 &3 


The e; are then called the principal dielectric constants. We 
now consider those solutions of Maxwell’s equations which 
can be represented by plane waves. Thus, all field inten- 
sities are to propagate as plane waves: 


E= Beth =—0. H= fe ke eitk: x—wt) 
wo @ 
k = i — 
WG Pe 
n= unit vector in the direction of propagation |, [12.9] 


v = phase velocity 
2 


c : ‘ 
p= rr n*, = index of refraction 


The vector n is also called the wave normal (since it is 
perpendicular to the surfaces of constant phase). We have 


Pp=n. 


We now substitute the forms given in [12.9] into Maxwell’s 
equations [12.1]. There follows 


cw to , cw ik= 2p, [12.10] 

D=—pxH, H=pxEk, [121] 
D = p*E— p(p-E) 

Ba nt ee ai) ‘ [12.12] 
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From Eq. [12.11] we conclude that D and H are perpen- 
dicular to p, and that D, E, and p lie in a plane. 
Equation [12.12] together with Eq. [12.3] yields 


il D 
— ae =. ; 
(- =) +n 13 2 ns) 0 [12.13] 
If we set 
f= Dy [12.14] 
k & 

then Eq. [12.13] becomes 

1 af 

n &; 


Equation [12.15] is equivalent to [12.13] if we require 
that A be taken so as to satisfy the requirement 


> DN = 0 .) [12.15a] 
Ee 


which follows from divD=0 and Eq. [12.10]. (This con- 
dition follows directly by multiplying Eq. [12.15] by n, 
and then summing.) Thus, instead of [12.13] we can use 
the system of Eqs. [12.15] and [12.15a]. 
We solve Eq. [12.15] for D,, 
An; Anp; 


re AE 12.16 
D; 1/n? —1/e; 1— p?/e; , [ ] 


and then take the scalar product with p. There follows 


Pi 
SS eS 0, UWA 
gee San, De aie [ J 


Written differently, 


_ oPi— Pp piles + P vile: _ 2 2 Piles 
2 te Bee, 2\ri +P 1— p*/e; 
pil &; 


pt PD a ale; pe, 
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Thus, 


Piles Bi 
Se ae 12.18 
2 peje —1 ars 
If we clear the denominator in Eq. [12.17], then there 
follows 


il aN ffl 1 1 ie 1 
2 ae ee Ee 2 aay capac ei. SS 
™ (<; A [- =) he Ms (- =| & | 


+n (=-2) (=< =0. [12.19] 


n= &)]\n* ~ & 
This yields a quadratic equation for 1/n?. Thus, in gen- 
eral, we obtain two different roots for the refractive index. 
To each wave normal then correspond two different phase 
velocities. 


Special case: &,=€; 

Form the preceding discussion (see the remark to 
Eq. [12.7]), this implies an optical rotational symmetry 
Equation [12.19] splits into two factors: 

iL 1 i 1 1 1 
(- = {mt (= sai + (ng +ni C= =) =0. [12.20] 


NW & nm? — Es) 1 

Case 1. The first factor is zero: 1/n?=1]e,. 

The corresponding index of refraction is independent of 
the wave direction n. This is the so-called ordinary ray. 
If we denote the plane through axes 2 and 3 as the prin- 
cipal section, then we can say that the D corresponding 
to the ordinary ray must lie in the principal section and 
be perpendicular to n. This is true because in Eq. [12.15] 
we must have A=0 (we set i1=2); for «=1, however, 
(1/n?—1/e,) 40, so that D,=0. Thus, for the ordinary 
ray E is parallel to D (because the latter is in the prin- 
cipal section). From Eq. [12.11], the magnetic vector for 
the ordinary ray then becomes 


H=n(nxB) = SenxD, [12.21] 
2 
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If n is perpendicular to the principal section (i.e., parallel 
to the 1 axis), then the direction of D in the principal 
plane (and thus the direction of the ordinary ray) is not 
defined. 


Case 2. The second factor in Eq. [12.20] is zero: 


ai(—=)+at+0n(5—2) =o. 
n= & 


nN & 
From this (since n?=1), there follows 


1st hte 


aa _ [12.22] 


If we denote the angle between nm and the 1 axis by 9, 
then there follows from Eq. [12.22] 
cos? sin? 1 
cai pk, [12.23] 
Eo &y n® 
Here, the refractive index is dependent on the direction 
of the normal. This is the extraordinary ray. 
We consider an extraordinary ray lying in the 1-2 plane 
(the 2 and 3 axes are suitably chosen), so that n,=0. 
Then, from Eq. [12.15] it follows that 


(= — =| D,;=0 and therefore D,=0. 
nm & 
Thus, the D vector lies in the 1-2 plane and is perpendic- 
ular to n. If we consider an analogous ordinary ray, then 
we know that its D vector lies in the principal section and 
is perpendicular to n. Thus, we conclude that the ordi- 
nary and extraordinary rays are polarized perpendicular 
to one another, With the use of [12.21], this can also be 
expressed by saying that H of the ordinary ray is in the 
same direction as D of the extraordinary ray. 

If in Eq. [12.23] we set y=0, then the ordinary ray 
can no longer be differentiated from the extraordinary ray. 
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Thus, in our special case there is one (and only one) di- 
rection for which Eq. [12.19] has a double root for 1/n?. 
Hence, crystals of this kind are called optically uniaxial, 
and the direction for which there is the double root is the 
optical axis. 


Geometrical interpretation of the general case. 
We set 
x = constant xD. . [12.24] 


From Eq. [12.16], then, 
“1 — 0 [12.25] 
If we take the ellipsoid 


> 1 [12.26] 


and intersect it with the wave plane [12.25], then the 
principal axes of the ellipse formed by this intersection 
are proportional to the squares of the refractive indices 
m, and n, which belong to those polarizations of the plane 
waves of wave~normal n whose vectors D fall along the 
principal axes under consideration. This statement is the 
Same as saying that the squares of the refractive indices 
m, and n, are proportional to the extrema of the quadratic 
form 


«i = extremum [12.27] 


+ 


under the requirement of the auxiliary conditions [12.25] 
and [12.26]. With suitably defined constants, the method 
of Lagrangian multipliers yields 


post yn =0. [12.28] 


If we take the scalar product of this with z,, then upon 
using Eq. [12.25] and [12.26] there follows for the extreme 
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value of [12.27] 


at 

Dae. ——. [12.29] 
i BX 

On the other hand, D satisfies Eq. [12.15] and is given 
by Kq. [12.24]. Substituting the latter as a geometric fact 


into the former yields the ellipsoid 
a th =0, [12.30] 


which is compatible with Eq. [12.28] only if ~=1/n?, so 
that 
(S24) a. = constant xn’. Q.E.D. 


13. RAY VARIABLES 


a. The ray index 


Instead of the wave variables used previously, we can 
also introduce ray variables. To define these, we proceed 
as follows. Equations [12.11] are relations between p, E, D, 
and H. For the Poynting vector we had 


Cc 1 
= — W=—(E- 2), F 
Z, Ext = D +|H|?) [13.1] 
We now define the ray velocity U by 


Cc 
=; ¢- [13.2] 


The second part of Eq. [13.2] simultaneously defines the 
ray index s and the ray vector t. (We have assumed |t|?=1.) 
For isotropic media, the ray index s is equal to the index 
of refraction n, and the ray vector t is equal to the wave 
normal vector n. For crystals, however, this is no longer 
the case. A relation between the ray index and the re- 
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fractive index is obtained as follows: 


S Be cExH _cExH 
W HE-D+|H)}  |H|? ’ 


c 
er 
s 


since for a plane wave the electric energy density is equal 
to the magnetic energy density. Indeed, from Eq. [12.11], 
there follows 


D-E=— (px H)- E= (px E)-H=|H|?. 
Thus, from Eq. [12.11] and [13.2] there follows 
c cExH cExH ct 


=" H- (px E) ~  p-(HxE) p:t’ 


from which p:t=s, or 
[13.3] 
b. Duality 


If, in addition to the above, we also introduce the 
variable 


1 
: | ae t ’ [13.4] 
then we can convince ourselves that all equations in this 
chapter remain correct if the following substitution of 
the first by the second line is made: 


Dpnr «s n H 


Eg? 1s" Ve, tH. (e) 


The variables in each column are said to be duals of one 
another. 

The substitution [13.5] produces a correct equation out 
of each equation, provided that the original equation is 
correct. Thus, we must show that Eq. [12.11] permits 
this substitution: 


c c C 
HxS= gg fi (Ex) = 7 E|H) — 7 H(E-#) =cWE, 
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and, analogously, Dx S=—cWH. However, with S=ecWg’, 
E=— q'xH, H—@q-D, [13.6] 


which corresponds exactly to Eq. [12.11]. From this, how- 
ever, everything else follows: from p we obtain q’ by sub- 
stituting n—>1/s, nt, and the equation (for principal 
axes) D,—=e,E, remains correct upon replacing E->D and 
&—>1/e,. Q.E.D. 

We know that the direction of the Poynting vector, and 
thus that of ¢ and q’, is the direction of the transport of 
energy. Thus, q’ is the direction in which a wave packet 
propagates, and n is the direction in which the phase 
propagates. Thus, all explanations of Section 12 can, with 
the aid of the scheme of [13.5], be directly taken over 
for the ray variables. 


c. Relation between wave normal and ray direction 


We have Eqs. [12.11] and their duals: 
i tx Pp, E=Hxg' . [13.7] 


From these, we obtain 


D-q'=(Hxp)-q'= (pxq’)-H, [13.8] 
E-p = (Hxq’):p= — (pxq’')'H, [13.9] 

so that 
p:E=—q''D. [13.10] 


On the other hand, we have Eq. [12.12], 
D = p’E — p(p'E) [13.11] 
and its dual 
E=q"D—q'‘(q''D). [13.12] 
In component form (for principal axes), 
p(p:E) = p’E;—«.F;, pene 
— qi(q':D) = E; _ qe; E, . [23.12'] 
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Dividing Eq. [13.11'] by Eq. [13.12] and using Eq. [13.10] 
yields 


Arrogant [13.13] 


Furthermore, we have Eq. [12.17], 


i ae 
»3 1— ple; me 
whose dual is 
re 
qi 


and Eq. [13.3], which can be transformed to 
>a 1 [13.15] 


From Eq. [13.13] there follows [we wish to calculate 
U t 
q,=9,(P)] 


(los Pi .- alee t Pi 
G= (a A) t= pg, = ing —1)-o i pe 
so that ~ 


pq: = p: + (p2q’?—1) [13.16] 


le ae 
or 


Ra, —— flere | 
Pd: — Di = (p*q ) i ae ae: 


When taking the vectorial square (that is, squaring the 
ith component and then summing) and using Eq. [13.15], 
this gives 


4 
= 3p" + pa (pig? — 1)" > Ps - 

G —p i) 
Thus, 


p? = (p*q” —1) ee Rea 
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If we set 


then from Eq. [13.17] there follows 
pq’ —1=p'/K, 


and this when substituted into Eq. [13.16] yields 


so that 


= ‘ |: [13.19] 


34 
i= pals, aT EK pie) 


We can, of course, without calculation write down the cor- 
responding dual equation (using substitution [13.5]): 


allel: il 1 
Pi = 4 Fe se a ’ [13.20] 


with 


At the same time, we see that 
hig? —= Kip. [13.21] 


Thus, Eq. [13.13] is solved for q, and p,. 
We can represent all of this in still another way by ex- 
plicitly introducing the function @ from Eq. [12.17]: 


== 5 ie a ae [13.22] 


Then, by differentiation we get 


oG 
; Op: = 1 ore 
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The second term in this is (see Eq. [13.18]) 


5 Piler sp {(L/ex—1/p*) + 1/p*} ve 


1— pes (1— p?/ex)? 
= = (-244 K)= Kip’. 
Pp 
Thus, we obtain 
oG Ds 
ee ee =a 13.23 
Ae Fie, gene [ ] 


and, by comparison with [13.19], 


<— = Ky}, [13.24] 


[13.25] 


(in accordance with [13.15]). Following the scheme of 
[13.5], the corresponding dual equations are 


tees os 
F(q') = a =0, [13.26] 
. oF 
———— ‘a. 13.5 
oF 
K'=Ya5n- [13.28] 


d. Geometrical interpretation 


1. In Section 12 we have already seen that the refractive 
indices corresponding to a given wave normal can be ob- 
tained from the intersection of an ellipsoid and a plane. 
In accordance with the scheme of [13.5], the ray indices 
can also be obtained in this way. If the ellipsoid } e,a?=1 


is intersected by the plane x-t=0, then the principal axes 
of the ellipse formed by the intersection are proportional 
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to the reciprocal of the two ray indices 1/s, and 1/s, that 
correspond to the ray direction ft. 

2. Equations [13.22] to [13.25] and their duals, Eqs. [13.26] 
to [13.28], also admit a geometrical interpretation. Here, 
the ray direction corresponding to a wave normal p is 
obtained by dropping a perpendicular from the origin to 
the tangent plane which passes through the point of inter- 
section of the wave normal p with the so-called wave sur- 
face G(p)=0. This is true because Eq. [13.24] states that 
the vector q’ is perpendicular to this tangent plane. The 
length R of the perpendicular is (see Eq. [13.15]) 


R= p:t=s(p-q')=s. [13.29] 


Thus, the length is exactly equal to the ray index. The 
pedal surface of the wave surface G=0 (i.e., the surface 


Pedal surface 


Figure 13.1 


consisting of the locus of the points of intersection of the 
perpendiculars with the planes tangent to the wave sur- 
face) is, therefore, the surface obtained by measuring off 
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from the origin a distance s, corresponding to the ray index, 
along every possible ray direction. Conversely, the surface 
G =0 is obtained as the envelope of the planes which are 
constructed perpendicular to radius vectors of length s 
(equal to the corresponding ray index). The dual construc- 
tions are, of course, also valid: the wave normal corre- 
sponding to a given ray direction is obtained by dropping 
a perpendicular from the origin to the tangent plane passing 
through the point of intersection of the ray vector q’ with 
the so-called ray surface F(q’)=0. The length Q of this 
perpendicular is 1/n, that is, exactly equal to the reciprocal] 
of the index of refraction. Consequently, the pedal surface 
of the ray surface F=0 (i.e., the locus of the points of 
intersection of the perpendiculars with the planes tangent 
to the ray surface) is the surface obtained by measuring 
off from the origin a distance 1/n, the reciprocal of the 
corresponding refractive index, along every possible wave 
normal direction. Conversely, the surface f=0 is obtained 
as the envelope of the planes which are constructed per- 
pendicular to radius vectors of length 1/n (n = correspond- 
ing index of refraction). 


3. There exists a relation between the vector q’ and the 
wave-kinematically defined group velocity: 
U=ceq'. [13.30] 
As proof, it must be shown that U satisfies the same equa- 
tion as q’. From our previous work, we have 


@ @ : @ k; 
k= nn =~ p= kop with i so that Die 


The group velocity U was defined by 


ow 
U= x @ = w(k,, ky, ky). 
From the above equations follows 
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For the proof, we introduce the variables k,, k,, k,, and k, 
(homogeneous coordinates) in place of the p; in the func- 
tion G: 

1 pi 1 ki 

G — ead ss eat = 

5 Da =a Di Reps, (=): [13.31] 
Since Eq. [13.31] defines a homogeneous function of the 
zeroth degree in k,,k,,k,,k;, then the Euler relation holds: 


ets et, [13.32] 


Now, k, must be such a function of the k; that for a vari- 
ation of the k, (transition to a neighboring wave normal) 
the form G must still vanish. This yields the condition 


OF Om (PF) _ 4 gp LeGF__ oh 
Ok, Ok; | \Okije, ky Op: Oty OF,” 


From this, upon using Euler’s relation [13.32], it follows 
further that 


+ 


is yg _ (BO ger, 20) 
Op: ° Ok Ok; =) Ok, \ko Op) 0k; — P | Ok; 
Comparing this with Eq. [13.24], we get 
Ok, 
qi Ok,’ 
or 
7 1eo Ut 
qi — e Ok; a eC U,; 0 Q.E.D. 


In concluding these paragraphs we may note that here the 
quantities p, q’, and G@ play roles analogous to those of 
the similarly designated quantities in geometrical optics (see 
Section 3b, pp. 16-20). 


14. SINGULARITIES 


There exist singular directions for which no unique ray 
direction is associated with the wave normal (and, of course, 
the same is true for the dual quantities). This is the case 
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when Eq. [12.17] or Eq. [13.31], which determines the 
index of refraction, has a double root. In that case (in the 
principal axis system) we have 


n 


96 2 ie ie, 


With the abbreviation 
il aL ff al i 
r-(-)G-26-2) 
n® &)/\n® &/\n® 6&3 
we have 


ae (=—- | ee 
ave = ni (= Ea) \N2 Ee. 1M nm? si) \nN® &, 


1. ieee 
2 
— = | 
TMs ( =) (3 a 
Now, if NG =0 is to have a double root, then it must also 


be true that (NG@)’=0, where the prime denotes differen- 
tiation with respect to 1/n?. This yields 


(definition of a double root) . 


NG@=0 
(N@)'=0 


Hence, 
NG — vp 


As can easily be shown, G’ cannot be zero. Therefore, 
for a double root in 1/n?, 


(0 Se 

nN? &)\nN® e/\n®  &, 
This is the case when, for example, n?=e,. Moreover, all 
the e; are to be different. It then follows that 


so that 
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If the preceding is substituted into the general formula 
for NG, then we obtain 


1 al 1 al 1 1 
2NG =|——— 4 2 
Fe S| (n (5 =) sia E =} 


For a double root, the quantity in curly brackets must 


also vanish: 
1 3 
a(d-d)on(t-2)-0 
€, =&3 Eg & 


This equation can be satisfied only if the two bracketed 
terms have opposite signs. Thus, 


&<&<€ (or vice versa), [14.1] 


and for the directions for which the double roots occur 
there follows 


~~ Vite l/é ~=/) Ee 


1/e, = 1/es E2(€3 — Ex) 


=i) 25% fe, _ s]) 

1fe,— Ley &2(€3 — ‘A 
Thus, there are two such intrinsically different directions 
which are called the binormals or optic axes. In accordance 
with the general construction, the refractive indices corre- 


sponding to a given wave normal are obtained by cutting 
the ellipsoid >'#?/e,=1 with the plane x-n=0. If n is 


[14.2] 


along a binormal] direction, then the intersection is a circle 
and the binormals pass through its center and through the 
umbilical points! of the ellipse in question. Since the 
circle has no principal axes, D and, therefore, E are inde- 


1 Translator’s Note: If we consider a family of parallel planes, each inter- 
secting a quadratic surface in a circle, then the points of tangency of these 
planes with the surface are called the umbilical points of the quadratic surface. 
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terminate. The only requirement is that D must be per- 
pendicular to the binormals. Now, from Eqs. [12.11] and 
their duals, we know that E, D, t, and n (in our case n = b) 
lie in one plane. (See Fig. 14.1.) However, D can have 


D 


Figure 14,1 


any direction perpendicular to b, so that the whole figure 
is undetermined with respect to a rotation about b. Thus, 
we see that € lies on a cone. It should be noted that the 
angle is dependent upon the azimuth, so that in general t 
does not lie on a circular cone. Thus, we have proved 
the hypothesis asserted at the beginning of this section: 
if the index of refraction is a double root of Eq. [12.17], 
then no unique ray direction is associated with a given 
wave normal. 

We now wish to derive the equation of this cone. In 
the principal axes system (which is arbitrary in a plane 
perpendicular to 6), D has the following components: 


D = {—Aby, pw, Ab}, [14.3] 


for then the scalar product with b is equal to zero. Thus, 
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E becomes 


E= |= 0,4, dele [14.4] 
€y &2 &3 


The generatrix of the desired cone has the direction of t 
and, therefore, is perpendicular to E and lies in the plane 
containing E and b. This leads to the requirements 


b-(Exx)=0, and E-x=0, 


where we have denoted the coordinates, in the principal axes 
system, of an arbitrary point on the cone by x= {£, H, Z}. 
In component form, 


b,(E,Z — E;H) = b,(E, H — E,=) =0, 
and 
E,S+ #H,H+#H;Z=0. 


Into these we substitute from Eq. [14.4]: 


and 


A 
mec 1 3.7 = ., 
&} Eg &3 


From these, we eliminate A/u: 


pe WDE 
(6,2 6,8)" =18(2+2), 


2 és Et 
He =; (> z— 22), 
Eg €1 €3 


and, by dividing, 


2 2 
(b,2 —b,5 (3- 22) == (242). [14.5] 
&j Eg €3 €; 


Equation [14.5] is clearly that of an elliptical cone. If 
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the b; from Eq. [14.2] are substituted, there follows 
= 5?(e, — &) + H*(es—&:) 


E3€1 


Equation [14.5] can be further transformed as follows, so 
that the geometrical form of the cone becomes clearer. 
From Eq. [14.5] there follows 
ae) es ‘ve ail 
te (2 + a)+e 2 +5ab),(->-=)+ 3:2 = F 


€3 


1 bs Ot 
(2+2) (S+H+ 2) — (2 2453) 


€3 & 


2 2 
(242) (Bt. Ht 4.20) — (2b, + 5b,) (22 +5%)=0. [146] 
3 1 
From Eq. [14.6] we see that the cone defined there has the 
same curve of intersection with the plane 


pee 0 [14.7] 


(that is, the wave plane) as does the sphere 


(2+ 2) e+m+2)—o(22+2%) <0 
€s &y €} &3 

with the wave plane [14.7]. We thus have the important 
theorem that the planes perpendicular to the binormals in- 
tersect the corresponding ray cone in circles. Furthermore, 
we know that the cone contains the vector b in its sur- 
face (we need only to choose D parallel to the 2 direction), 
and that the entire geometry is symmetric with respect to 
the 1-3 plane. The angle «’ between b and the other gen- 
eratrix lying in the 1-3 (symmetry) plane, which is a meas- 
ure of the opening of the cone (clearly, we can not speak 
of an opening angle for the case of an elliptical cone), is 
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also the angle between E and D in the 1-3 plane. How- 
ever, this latter angle can be easily calculated: 
tan «'= tan[(E, D) (in the principal plane)] 

= + tan[(E, 1) — (D, 1)] 

= [tan (E, 1) — tan(D, 1)]/[1-+ tan (E, 1) tan(D,1)], 

tamn(D; 1) =D /D, = — b,/b; 

(since D is perpendicular to b), 

tan (E, 1) E,/E, = D,e,/(Dyé3) = — 0,81] (ses) , 


so that 
— &b,b,-+ €5b,b, 
£05 = €,D} , 


+ tan «’=-+ tan (E, D) = 


If we substitute 6, and b, from Eq. [14.2] and order the 
terms, there follows 


tana’ = ie pee 21) [14.8] 


E1,€3 


The entire discussion can be given in a completely analo- 
gous way for the dual quantities. The most important 
theorems are then as follows. 

There exist singular directions for which there is no 
unique wave norma] associated with the ray direction. 
This is the case when the equation that determines the 
ray index has a double root. The ray directions for which 
this is true are called biradials or ray axes, and their di- 
rection cosines are 


— E3— &2 
n= a ae’ r,=0, naa Or “a 
et 3— €1 


The associated wave normals lie on a cone which inter- 
sects the planes perpendicular to the biradials in circles. 
The angle f’ (a measure of the opening of the wave cone), 
which corresponds to the angle «’ discussed above, is 
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given by 


ja /@ — Sa)les = 64) | [14.8] 


Eo 


These results can also be obtained without calculation from 
the scheme of Eq. [13.5] and the corresponding formulas 
for the binormals. 


15. LIGHT ENTERING AND LEAVING CRYSTALS 


Consider the case shown in Fig. 15.1. Here too, we have 
the law of refraction: 


nN, Sin «, = nz, Sin «, 


Ne = Nq(%) Ba! 


The system of equations [15.1] is very complicated to solve 
analytically. Huygens, however, gave a construction with 


m,: isotropic medium 


Ns: crystal 


Xs 


Figure 15.1 


which it is easy to get an overall view of the various re- 
lationships. Let A,B, be an edge view of part of a sur- 
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face of constant phase of a plane wave incident on the 
crystal (the surface of constant phase is perpendicular to 
the plane of Fig. 15.2). The construction is such that the 


: erystal 


n,: isotropic 


medium 


Figure 15.2 


optical path from B, to B, is the same as that from A, 
to A,; that is, the surface 2 is so defined that the length 
of each line drawn from B, to 2 is the reciprocal of the 
index of refraction associated with the direction of that 
line. As can easily be seen, the plane perpendicular to the 
plane of the drawing which passes through A, and is tan- 
gent to the surface 2 is a surface of constant phase (the 
optical path from C, to C, is always the same, no matter 
where C, is chosen on A,B,). Of course, the point of con- 
tact does not necessarily lie in the plane of the figure (as 
it is drawn here for simplicity). Since the wave normal is 
always perpendicular to the surface of constant phase, it 
is obtained by dropping a perpendicular, Ben, to the 
tangent plane mentioned above. The associated ray direc- 
tion can then be found with the use of a previous con- 
struction. It should be noted that the surface 2 consists 
of two sheets; thus, in the crystal there are two wave 
normals and two ray directions for one plane wave incident 
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on the crystal. The plane wave is split into two, as is 
the light ray.? 

A singularity occurs if the normal direction (B,B,) cor- 
responds to a binormal of the crystal. This results in 7n- 
ternal conical refraction. The refracted rays then lie on a 
cone (see Section 14). If a light ray is sent through a 
crystal with two plane parallel surfaces at such an angle 
that internal conical refraction occurs, then the ray emerges 
as an elliptical cylinder. 

The corresponding dual phenomenon is external conical 
refraction. If a light ray which is parallel to a biradiai 
direction within a crystal emerges from the crystal, then 
the emerging wave normals lie on a cone produced by the 
refraction of the indeterminate wave normals which are 
associated with the biradial ray and which lie on a cone. 
A suitably chosen convergent pencil of light incident on 
a crystal plane will, as a result of this phenomenon, form 
@ single ray within the crystal and will emerge from the 
other side as a divergent pencil [A-5]. 


* Translator’s Note: This phenomenon is known as birefringence or double 
refraction. 


Chapter 5. Molecular Optics 


16. DISPERSION BY UNDAMPED OSCILLATORS 


We now make specific assumptions about the medium 
in which the optical phenomena take place. In this manner 
we shall attempt to explain the optical constants. 

We assume that charged particles are distributed through- 
out the medium under consideration. Usually it will suf- 
fice to replace the discrete distribution by a continuous 


one. If 
a=, [16.1] 
k 


in which > is the sum over all particles in a unit volume, 
k 
and ¢, is the charge of the kth particle, then we denote 


P= Dd eX [16.2] 
k 


as the electric dipole moment per unit volume. The cur- 
rent density is 


i=P=>e,0,. [16.3] 
k 


For many purposes it may be assumed that the positive 
and negative particles are of the same type, so that we 
can write 

P= N4.64.%4.-+ Nee2 x. ; [16.4] 


where W,, ¢,, x, represent the number of particles per 
unit volume, the charge, and the position of the center 
of charge, respectively, of the two kinds of particles. 
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For this kind of a medium the Maxwell equations are 


curl E=—=H, div H=0 
; [16.5] 
Ciel =D: divD = 0 


To these must be added the relation between D and E: 
D=E-+4aP. [16.6] 


Here, as before, we have assumed “4=1 (no magnetization). 
For a wave field periodic with the time factor e"*”, we 
have already seen that D=7n’E, so that with Eq. [16.6] 
there follows 

4nP = (n?—-1)E. [16.7] 


Then, the Maxwell equations have as solutions waves with 
phase velocity c/n. 

Furthermore, we will assume that the charged particles 
of the medium behave as quasi-elastic oscillators with 
respect to external influences. Then, the restoring force 
on a particle is 

2 K=—fzx, 


where x is the displacement from the equilibrium position. 
Without external forces, the differential equation of mo- 
tion for the particle is 


m*e+fx=0, ma, =f. [16.8] 


This assumption is consistent with the occurrence of sharp 
spectral lines. 

Now, assume that an electric field acts on the medium. 
This gives rise to an additional force 


K, =e [16.9] 
on a particle, so that 


mx + wex) = eE. [16.10] 


Sec. 16 | DISPERSION BY UNDAMPED. OSCILLATORS iil 


If the electric field is periodic, that is, if 

Ee, [16.11] 
then Eq. [16.10] admits a solution of the form 

B= 46-18 , [16.12] 


Substituting Eq. [16.12] into Eq. [16.10] yields 


m(— w?-+ w¢) x, = eE, , [16.12"] 
so that 
e 1 
and 
i 
= E. [16.14] 


m (ws — w) 


From this, we obtain a formula for the electric polarization 
and the index of refraction. If N, is the number of par- 
ticles per unit volume having a free oscillation frequency ap, 
then there follows (if there is no polarization in the ab- 
sence of the electric field E) 


Pay og 16.15 
aa om (we = w*) 9 [ es ] 

and from Eq. [16.7] 
re | ee a ees [16.16] 


m (we — w*)" 


If there are several types of oscillators, say N, of fre- 
quency «, per unit volume, then there follows 


e2 1 


reat [16.167] 


k 
Here, of course, we have made the essential assumption 
that the number of particles within a cube the length of 
whose side is a wavelength is large. This is no longer cor- 
rect for the case of x rays. 
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Equation [16.16’] is only conditionally correct because 
the field intensity E is not equal to the field intensity E’ 
that produces the polarization. The field intensity in 
Eq. [16.10] is understood to be the exciting field, while 
that in Eq. [16.7] is understood to be the average field. 
The connection between E and E’ is determined in electro- 
dynamics in the derivation of the Clausius-Mosotti equa- 
tion. We have 


E = E'—9gP, [16.17] 
where for isotropic bodies 
gy = 42/3. [16.17'] 
With this correction, there follows 
4nP = (n?—1) E’'— - (n?—1)P 


or p | . [16.18] 
= (n? +2) P = (n?—1)E' 


This, now, is the equation which must be substituted into 
Eq. [16.15]. There then follows 


vl e 1 
= 16.19 
Ww+2  °m(o—o)’ [ ] 
or, if we again allow several kinds of particles, 
ni —1 ek al 
3 == 4 N,— 16.19’ 
m+ 2 m2 * my (wt — w*) oe || 


Of course, even this equation is not correct since the 
right-hand side of Eq. [16.19] has a pole at w=a,. Thus, 
we will have to assume a damping force later. For the 
present, we wish to pursue somewhat further the results 
that follow from this model of quasi-elastic, charged par- 
ticles. 
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a. Molar refractivity 


1. Let the moving particles be electrons, let M be the 
molecular weight of the substance under consideration, and 
let @ be the density. We define ‘‘molar refractivity” R by 


= nt 2 0 %s [ Gg J 
Then, there follows from Eq. [16.19], 
4a ée ie 
k=—A>— —_.. ; 
3 kM, (wi = w”) L16 pa 


Here, A is Avogadro’s number, and f, the number of oscil- 
lators per actual molecule. According to our model, it 
would be expected that the /f, are integers (since they 
are considered as being particles that oscillate). In real- 
ity, even for the hydrogen atom (a single electron) there 
are many such “‘oscillators,’’ and only the sum of all these 
‘“‘numbers”’ of oscillators is equal to the number of elec- 
trons Z in the molecule: 
>i =Z.- [16.22] 
k 
Thus, f, is called the oscillator strength corresponding to 
the frequency w, [A-6]. 
For the limiting case w> w, for all k, there follows from 
Eqs. [16.21] and [16.22] 


e 
mo? 


4n 


(Here, ¢ and m are constants of the electron.) In this 
limiting case, the electrons are considered as being free. 

2. Let the moving particles be atomic nuclei. For oscil- 
lations of diatomic molecules we have f=1 and 1/m 
=1/M,+1/M,, where M, and M, are the masses of the 
nuclei. Since the oscillator strength is an integer, the model 
of oscillating particles does not break down as it did in 1. 
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b. Conservation of energy 


We again neglect the difference between E' and E. The 
Poynting vector is 


S=—3ReH, 
An 
so that 
div S = 7, (Hr eurl E— E-curlH). 


However, since we have the familiar energy equation 
divS+W=0, 


then if we use Maxwell’s equations [16.5] and [16.6] there 
follows 


W= (HY E-E) + E-P. 


However, in a unit volume we have 
P= > Weer Xx 
k 
and 


he e€ 
x, + wo; xX, a = E ; 
so that there follows 


1) > Wm, Xe° (He-+ wh *,) 
k 
d Mars 
= ay DMF t+ ata. 
Thus, after integration, 
W=— (H+ Es) 4 5, (a2 4 what 
= aR sil en by (Fe + Exe), 


Wrrar aa Ue mag. a Wag 


terial oscillators ? 


which is the law of conservation of energy. 


Sec. 16 | | DISPERSION BY UNDAMPED OSCILLATORS 115 


Here too we can show that 


is the group velocity. In order to see this, we first cal- 
culate the energy density of the material oscillators. We 
have just shown that this is 


™m ° 
Wrst a 20a (ag + wz XZ) g 
k 


Now, however, we have Eq. [16.14], so that 


= (€/m)x = (e/m);, —iwt—id ’ 
ae aR re yeast Re E,e . [16.14'] 


Since we are dealing with quadratic expressions, we must 
first take the real parts of the complex quantities. When 
Eq. [16.14'] is substituted into the relation for W_,,, there 
follows 

m, _(e/m)e 

2 (wy — w*)* 


x {(Re (— iwE,e-'-**))? + wy (Re E, e~#-#)?} 


Wrst — > N, 


x {w? EG sin? (wt + 4) + wz EG cos” (wt + 6)}. 


For the time average we have the familiar result 


|E|? =|, |, 


from which there follows 


Thus, if we recall that |H|?=7?|E|?, then there follows for 
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the total energy density 
1—|. ee ato? 
Woe = =|)? (14m dor 3 are 


According to Eq. [16.16’], 


ee ii 
71S ee 16.16'5 
w—1 A eae [ J 
so that the sum in the last expression is 
Ce Oe +a" yr Ge 200? 
4n > N,— —, = n?—1+42 > VN. —- ———saa 
* 2M, (at — oF} +42 2M, hay 
d 
=n—1 —— (n?—1 
ni —1+ oa | ), 


which when substituted into the relation for W yields 
—— 1—— d 
Se ee 2 2 = (ro S 
W 3,1 = {2" + Oct y 
The Poynting vector S is 


c 
- == ee | 
S re | x | 
so that 
C 
|S|= Te n|E|?, 
and, furthermore, 
Wa 
|S| an dw 
On the other hand, we have U=0w/0dk for the group ve- 
locity. In an isotropic medium U=dw/dk, with w—=ck/n. 
Thus, 
ee dk da 
U do da 2) 
from which 


U =|S|/W = do/dk . Q.E.D. 
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17. DISPERSION BY DAMPED OSCILLATORS 


We start from the original form of the dispersion law 


1 


2 
a 
mM (wo — w*) 


[16.16] 
This gives us a singularity at o=,. Since this, of course, 
cannot physically be the case, we must introduce a damp- 
ing force so that the differential equation [16.10] now reads 


mx+fxtgx=ecE, ivi.t | 
which, with 
f=mo, and g=my, 
becomes 


B+ obx-+ye =—E. [17.2] 


We make the same exponential hypothesis as in Eqs. [16.11] 
and [16.12] so that there follows 


x(w8 — w — iwy) = <E [17.3] 


This, however, is exactly the same equation as [16.12’] 
but with w? replaced by w?+iwy. Thus, all conclusions 
remain the same if this substitution is made everywhere. 
Then, for example, Eq. [16.21] goes over into 


An é 1 
3 ee l= Sap ie 
Thus, there is no longer a pole at w=w,. If the index 
of refraction is computed from Eq. [16.20] using [17.4], 
then it is seen that the index is complex. However, in 
accordance with the previous work, this indicates an ab- 
sorption. From expression [17.4], it is seen that R? has 
an extremum (for each oscillator) when the denominator 
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becomes minimal: 
| w8 — w* — iwy|* = (ws — @”)? + w*y” 


which is to be minimal for #=o,,. Then, 


m* 


— 2 (02 ~ 02) 20m + 20my* = 90, 
so that 
wor = we—y*/2. [17.5] 


If there is no external field E acting on the oscillator, then 
it will perform free but damped oscillations. For describ- 
ing these oscillations we set 


& = aoe Rutt [17.6] 


If Eq. [17.6] is substituted into the differential equation 
[17.2] with E=0, there follows 


(— iw'— pw)? +5 + y(—iw’—p) =0. 


We separate this into real and imaginary parts. The im- 
aginary part is 
2iw'p—iwo'y =0, 
so that s 
w= yd. | [17.7] 


The real part (with Eq. [17.7] substituted in) is 


—w"+y"/4+ wo — y?/2 07 
from which 
wo =op—y"f4. °° | [17.8] 


Thus, the frequency of the free damped oscillation is shifted 
somewhat with respect to the eigenfrequency w, of the 
oscillator. 

From the theory of linear differential equations, we ob- 
tain the most general solution of Eq. [17.2] by the super- 
position of a particular solution, Eq. [17.3], and the gen- 
eral solution, Eq. [17.6], of the homogeneous equation. 
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The damping has the consequence that the oscillations 
do not become infinite at resonance. 

We obtain the physical significance of the constant y by 
considering the energy of the oscillation. As is well known, 
the energy F of an oscillator is 

m f 


eg Siar) aa 
Es ee ae + wpa?) . 


The time derivative of this energy is 
E = mix(% + wie) . 


If we assume free oscillations, then from Eq. [17.2] there 
follows 
E=—ymz?. 


However, for the time average (over one period) we have 


so that 
we he [17.9] 


which gives the significance of y. 


Application to the optics of metals 

We use a model of a metal in which the electrons move 
freely in a rigid ion lattice. Thus, we must set f=0, g40 
in the fundamental differential equation [17.1] which then 
becomes 


mx +g9x=ecE, [17.10] 
= ius eee 

=—E. We Ika 
K+ yx me [17.10’] 


In electrodynamics, the case where E is constant in time 
is considered. As the stationary solution (x= 0), one ob- 
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tains 


. € 
x=v-——E 


my 
i— Ney 0,04. pea) 
Ne? 
op = 
my 


In optics, ¥ must, of course, be considered too. The so- 
lution of Eq. [17.10'] is obtained in the following way. 
We set w,=0 in all previous equations, and neglect the 
corrections between E and E’. Then, if in addition we 
replace w? by w?+ iwy following the general scheme of this 
paragraph, there follows from Eq. [16.16’] 


al 


2 
a 
m (— w? — wy) 


[17.12] 
The index of refraction n’ is complex, as is required by 
the optics of metals. For metals the optical constants 
were defined as follows (see Eq. [10.2]): 


< 4no 


no” =e+i—. 
Ww 


By comparing this with Eq. [17.12], there follows 


ail = 4a eg z 
m (w?-+ y?) 
wea OY 
o=N as) ; [17.13] 
°° (at +9") 


The quantity y may be frequency dependent. It can be 
determined theoretically only with a quantum-mechanical 
model. 


For free electrons the damping is, in part, produced by 
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collisions with fixed ions; this damping decreases with de- 
creasing density. However, even for ideal gases there is 
the so-called radiation damping, which is directly related 
to the fact that an oscillator radiates energy. 


18. SCATTERING OF LIGHT 


a. Fundamentals from electrodynamics ) 


We start from Maxwell’s equations, Eqs. [16.5], and the 
additional relation for polarization, Eq. [16.6]. We con- 
sider all terms to be oscillatory with the time factor e*” 
(that is, 0/Ot~ —iw), and thereby obtain (for vacuum) 


eurl E =7 = a 
: [18.1] 
4 
curl H+ io E aa 
In terms of scalar and vector potentials we have 
.@ 
H— curlA, E=i-A— grad ® 
[18.2] 


divA —i—=0 
From these, there follow the second-order equations 


w? 
V2@ + hd = 4n div P 
mn [18.3] 


The corresponding formulas for a medium with refractive 
index n are obtained if ¢ is replaced by c/n and H by H/n 
(then Maxwell’s equations remain correct). Then, from 


1See W. Pavuti, Lectures in Physics: Electrodynamics (M.I.T. Press, Cam- 
bridge, Mass., 1972). 


122 MOLECULAR OPTICS Chap. 5 


Kq. [18.1] we have 


curl E = i= H 
: [18.1'] 
curl H+ ““ nE = — = ion? P 


The vector potential A should be A/n in the medium. 
That is, we demand that H=curlA hold in the medium 
as well. This, however, is arbitrary; it could be done dif- 
ferently. Thus, from Eq. [18.2], 


= curl, = i= A— grad 


ae [18.2’] 
GivA— "= n2@ — 0 


If we introduce the abbreviation k=wn/c, then we can 
also write Eq. [18.2’] as 


H=cuwlA, E =i A—graad 


[18.2”] 
div A — ikn® = 0 
Furthermore, Eq. [18.3] becomes 
V’D + 2D = 4n div P 
Vi A ect dn P|? [18.3] 
where we have set 
Pe [18.4] 


If wo introduce the Hertz vector Z, then we have 


A=—iknZ, @=—divZ. [18.5] 
Furthermore, 
0p = — divP [18.6] 
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(polarization charge density), and 
V’Z+hZ=—4nrP. [18.7] 


The field intensities expressed in terms of Z are 


H=—ikneurlZ 
E=kZ-+ grad divZ ‘ ees) 
The solution of the differential equation [18.7] is 
etkr 
vhs =|P. : ae [18.9] 


In the special case that the wavelength /A is small com- 
pared with the linear dimensions of the region in which P 
is appreciably different from zero, the familiar Hertz oscil- 
lator is obtained: ? 


R~|xp— xQ| 


=| Pas, TT, = Il | [18.10] 
eErkk 
Z,=i R 


where P is the field point, Q the source point, and dV a 
volume element. With 


_ (%#p— Xo) 
Larceny ew 
where 0 is a suitably chosen point in the source region, 


we have 


ik ikR 
H= n G +5) gem) 3 
ee wk 1 2 ts.1t) 
E="7, {Hr +5 a 


3ik 3 
+ n(n-T) (-% -F+3| 


2 See PauLi, Electrodynamics. 


124 MOLECULAR OPTICS | Chap. 5 
The rate of energy radiation by a Hertz oscillator is 


dH 2 12 on (@\ 
— a =o | =5en() i § [18.12] 


b. Lorentz’s theory 


In his theory of scattering, Lorentz considered a homo- 
geneous medium. In accordance with statistical mechanics, 
density fluctuations occur which cause scattering. The en- 
tire medium is divided into volume elements which are 
small compared with a cube whose edge is a wavelength 
long, but which, nevertheless, contain many molecules of 
the medium. Because of this assumption, this theory fails 
for x rays. 

If N, is the number of molecules in a volume element V 
and A is the wavelength, then 


Sosa and Via. 


Statistically speaking, we have 
AN, AN, =0 for VV’. 


Consequently, the light intensities originating from dif- 
ferent volume elements may be added. 

If we denote an average refractive index by 7, then when 
considering the fluctuations Eq. [18.1’] becomes 


curl H = {n? + A(n?)} (- = E) ; [18.13] 


where we have assumed that, on the average, the medium 
is not polarized. However, the part of E due to A(n?) can 
be written as a polarization; indeed, upon comparison with 
Eq. [18.1’] this is 


1 A(n?) il 
a E, r= 7 AME. [18.14] 


PS =r 


Then, however, in accordance with the formula of Eq.[18.12], 
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the volume V radiates energy: 


dk 2en/w\' 1 ——_" 
— EF (=) aaa ENO. 


This is because II= PV. If the expression for the inten- 
sity, that is, 


c 
I — an n | E|? 3 
is introduced into the last formula. then there follows 


I (An). [18.15] 


We assume that the refractive index of the scattering 
medium is only slightly different from 1: then we can 
neglect the correction between E’ and E and can use 
Eq. [16.16]. If we collect the nonfluctuating terms into 
the constant «, then we have 


n?—1=aN,/V, [18.16] 
A(n?) = a 7 
(A(n))? = = (AN,)?, [18.17] 


and, substituting into Eq. [18.15], 


dE 1o _—~— 
—— =I— — a (AN,). 18. 
dt 6x cf a7) [igs] 
Equation [18.18] gives us the energy scattered per unit 
time from the volume V. If we consider a light ray. then 
its intensity loss along an interval dz must certainly be 
proportional to the energy scattered in that interval. 


Therefore, 
=— —V [18.18] 
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as can easily be seen. With this, from Eq. [18.18], 
4 [~~ 
—— = -— ——.4?— (AN,)?. [18.19] 
4 u 


It is customary to introduce the scattering coefficient s 
and to substitute the expression for « given by [18.16] 
into Eq. [18.19]. From this there follows 


aa ar, [18.20] 
with 


(olga 


5 = 2 a1) 1 1 (AY,) 


where we have used the abbreviation N=N,/V. 


From statistical mechanics (the theory of fluctuations) 
we have the formula® 


= kT 
OM ape, Avi 
or 
(ANy)? _ Nk? 
Ny VER: 


If this is substituted into the equation for the scattering 
coefficient, there follows 


= SO Kae 
3 A \= VPM y, « [18.22] 
For ideal gases 
Op 
ai ee = 
(sr), <= er 
so that for these the scattering coefficient is 
BEY = 
=n Wa’ [18.23] 


Formula [18.23] permits Avogadro’s number to be deter- 


*See W. PAULI, Lectures in Physics: Statistical Mechanics (M.I.T. Press, 
Cambridge, Mass., 1972). 


Sec. 18 | SCATTERING OF LIGHT 127 


mined. Furthermore, we see that short waves are scat- 
tered more strongly than long waves; the blue color of 
the sky is a result of this. 


ce. Scattering of x rays 


We assume that the wavelength / is small compared with 
the dimensions of an atom. Then, we have incoherent 
scattering from the electrons. Let S be the energy scat- 
tered per unit time, per unit volume. Then we have (radi- 
ation from moving charges‘) 
2 e 


c= 


a wNZ, [18.24] 


in which WN represents the number of atoms per unit vol- 
ume, Z the number of electrons per atom, and v the ve- 
locity of the electrons. 

Now, the frequency of light is much greater than the 
eigenfrequencies of the electrons in the atom. Thus, we 
set wm) (and w,) equal to zero. This results in assuming 
(almost) free electrons. For such electrons, 


mo = cE. [18.25] 


We substitute Eq. [18.25] and the relation for the incident 
intensity, I = (e/42) E*, into Eq. [18.24]: 
8x e! 


a 3 mict 


INZ. [18.26] 


The scattering coefficient is given by s=S/I (this follows 
simply from the definition in Eq. [18.20] together with 
[18.18’]). Therefore, 


2 \2 
= (=) NZ. [18.27] 
3 \me? 


Now, N can be expressed in terms of the density o, the 


4 See Paunl, Electrodynamics. 
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molecular weight M, and Avogadro’s number A: 
N=o0A/M. 


If this is substituted into Eq. [18.27] and if one then di- 
vides by the density, then the mass absorption coefficient o 
is obtained: 


Sx { e \2 A Z 
— —Z=04—. 18.28 
le) 284g 88 
Equation [18.28] was established by Thomson. The num- 
ber of electrons in an atom Z can be determined with it. 
The dispersion formula for x rays is obtained by set- 
ting w, equal to zero in Eq. [16.16]: 


e2 


n?—1=—47NZ . 
mo? 


This shows the possibility that x rays may be totally 
reflected. 


Refinement. We consider a polarized wave, whose E vec- 
tor is parallel to the z axis, incident on an atom. As be- 
fore, the electron is considered to be ‘‘almost free’; how- 
ever, we now want to consider the fact that not all of 
the Z electrons of the atom that scatter the radiation are 
located at O, the position of the nucleus. The effect of 
the incident wave is to produce a periodic acceleration of 
the electron, in accordance with Eq. [18.25]. However, 
when accelerated periodically, the electron produces an out- 
going (scattered) wave whose field intensity E, at a dis- 
tance r (from the electron) is calculated in electromagnetic 
theory: 


E,=—'— 0, ° [18.29] 
The notation | is to be understood as denoting the com- 


ponent of the acceleration v which is perpendicular to the 
direction of propagation under consideration. Consequently, 
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the magnitude of the field intensity E, at the point P 
(a distance y from the electron, making an angle 6 with 
@ axis) is 
e? etkr 
= mor ee) a cosé. 

Here, E(x) is the magnitude of the incident field inten- 
sity at the point x and E,, perpendicular to the direction 
of scattering, lies in the plane containing this direction 
and the direction of oscillation of the electron. If we as- 
sume that r> |x|, then the distance R= \PO| becomes 


R=r-+(n-x)+..., 


where n denotes a unit vector in the direction of the scat- 
tering. Hence, we have (as an approximation) 


e? ikR 


Ep= —'E,(x) eis o cos6 . 


If we write the wave vector of the incident wave as 
ky, 
then, in our approximation, 
E,(x)= E,e**, 
where E* denotes the field intensity of the incident wave 
at the origin. From this there follows 


Go es ke Cee 
a mot eee (ke—k) x —R . [18.30] 
If we assume that there are several electrons present in 
an atom, then we must write 


Es= 0080, [18.30'] 


(Ge 
"BS eithke—k) xn 
at 2 R 


where the sum is taken over the Z electrons. Of course, 
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—iwt, 


E® is periodic in time with the factor ¢ 
Fo = Fe-it+d) | 


The intensity of the incident wave is then 


ae oe JP % 
_ 4n 
The intensity scattered from the atom into the solid an- 
gle dQ is 


2 \2 
iT (=) cos?6|K|#dQ [18.31] 
me 


The factor 
KC= >: et(he—k)-X_ [18.32] 


is called the structure factor. 

If we have light which is polarized at an angle with 
respect to the plane defined by the incident and the scat- 
tered directions, then we must decompose its electric field 
vector into components parallel and perpendicular to this 
plane. The component parallel to the plane is scattered 
according to Eqs. [18.31] and [18.32]; the perpendicular 
component is scattered according to the same formulas 
evaluated at 0=0. From this we see that the entire cal- 
culation for this component is the same as before, except 
for the single difference that in Eq. [18.29], &, is the same 
as ». That is, cos@ is to be replaced by unity. 

The square of the magnitude of the structure factor can 
be calculated: 


[KP = KOK = 5 elthe bee g-ith. bem 
nym 


=Z+4+2 > cos(k(n,—n)-(x,—%m)], [18.33] 


n<m 


where Z represents the number of electrons in the atom. 
Now, if 


k(n,— n)+(X.— Xm) > 1, 
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or 


(,— M)*(X_z— X_) > a ; [18.34] 


then all possible phases occur as arguments of the cosine 
and, consequently, the sum is small. Thus, 


|K|?~Z. [18.35] 


If the incident light is unpolarized, then all directions of 
E? are equally probable and the incident intensity is di- 
vided equally among the two states of polarization. Thus, 
Eq. [18.31] can be integrated. There follows 


| (Gea) contoz ans [ (= :) aa 


1 e2 2 
=a | Jcost# sina dé + [sino a6) 20, 
2 me 


82 [{ e? 
Z 
I=> (aa) ie 
and this is the previous formula of Eq. [18.27]. The in- 


equality [18.34] is always unsatisfied within a certain cone 
which, however, becomes smaller as A decreases. 


so that 


Remarks. According to the classical theory as well as 
to the quantum theory, the electrons in an atom are not 
fixed but, instead, have a certain probability of being at 
a given position. For a particular electron configuration, 
let this probability be 


Wiwinen resided Vy, 


[--[war, ge 


Then, the structure factor is to be taken as a statistical 


with 
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average: 
xr=f... wiz hires ma, ... AVe. [18.36] 
If W separates into a product, 


W (x1, 02+) Xz) = Wy(X1)We( Xe) ... We( Xz) 


fossyav— ay 


then expression [18.36] can be simplified. If we define the 
abbreviation 


with [18.37] 


(ee | Wal Xpethe bend Vy, [18.38] 
then 
KP= > DPPL+Z=2+ 0 (FFL + EF). [18.39] 
nom n<m 


If, in addition, w, depends only upon r,, 


Wy{ Xn) = WalTn) , 


that is, if w, is spherically symmetric, then Eq. [18.38] 
can be evaluated further. With 6 representing the scat- 
tering angle, 

[&,— k|= 2k sin 6/2 , 

“a= > i k), x} ’ 

cosa=—u, 
so that: 
dV=r'dr2x du , 


B= [ar yee sin (2kr sin 8/2) 


2 y? 
con 2a ar . [18.387] 


Thus, for spherically symmetric distributions, F, is only 
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a function of (47/A) sin 6/2: 


4 
(eae (F sin/2) [18.38"] 


Here, if A again approaches zero (more exactly, (sin6/2)/A 


->oo), then F, approaches zero and, as a result, |\KP 
approaches Z, which is the same result as before. 

In the above considerations we neglected the fact that 
(according to classical theory) a Doppler effect occurs dur- 
ing scattering. This effect shifts the frequency of the scat- 
tered radiation somewhat with respect to that of the in- 
cident radiation. 


19. OPTICAL ACTIVITY 


a. Phenomenological theory 


The optical activity of a substance consists of the rota- 
tion of the plane of oscillation of polarized light upon pass- 
ing through the substance. This effect can be described 
by saying that the refractive indices for waves of opposite 
circular polarization are different. (The superposition of 
two waves of opposite circular polarization of equal am- 
plitude and frequency clearly results in a linearly polarized 
wave.) The angle y through which the plane of oscillation 
is rotated is obviously 


Eo Steet). [19.1] 


Le 2S 


where I is the path traveled in the optically active sub- 
stance, and n, and n_ are the refractive indices for right- 
hand and left-hand circularly polarized light. 

In order for optical activity to occur, it is of course ne- 
cessary that the optical equations no longer be invariant 
under reflection. This results from the fact that certain 
molecules are not mirror-symmetric (for example, sugar 
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solution, tourmaline crystal), so that the eigenfrequencies 
for right-hand and left-hand circularly polarized light are 
different. 

Ordinary (polar) vectors are invariant under reflections, 
whereas axial vectors, such as the vector products of 
ordinary vectors and the magnetic field vector, are not. 
As shown in Electrodynamics,® these vectors are actually 
skew-symmetric tensors which change sign upon reflection. 

We take the following (Maxwell) equations from electro- 
dynamics: 


curl E = ==5 | 
: [19.2] 
nb eto fe | 
c c 
H = B—421M 
I, = P+ ccurl M | ° Le 
Furthermore, we assume 
os 
= E— «a, curlE 
An 
—_— Bay ; [19.4] 
aa c 
a+a,=a 
There then follows 
pact ae ‘ 
L, = E E—acurlE 
4m ; [19.5] 
D =cE—4nacurlE 
where D satisfies the relation 
yee 
curl B = me : [19.6] 


With these equations, the following results obtain for a 


“ibid 
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plane wave (V~ik, d0/ot~—iw): 


E— E, e** *z—wt) ; 


with 
k= 2 nn. 
Cc 
With the abbreviation 
= dan =e [19.7] 
there follows 
D=cE+ifnxE. [19.8] 


Now, as is well known, the Maxwell equations yield 
1 oe 
curl curl E = —— D 
G2 


and so, for a plane wave (where E-n= 0), there follows 
—vnx(nxE) = —n?[(n-E)n—E]=D, 
so that with Eq. [19.8], 
WE =cE+ifnxE. [19.9] 


We now locate n along the z axis (n= {0,0,1}) and E 
perpendicular to it (E = {E,, E,, 0}). Equation [19.9] yields 
nE, = cE, — ifE, 


mE, = cE, + iE, en 


In the usual way, we introduce H, and E,, in place of 
iE, and E,: 


Sm sones a: [19.11] 
Ey = H,—1tHy, ; 
Then, 
te ae [19.10] 


{n? —e—f} En =0 ie 
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or, written explicitly, 


(nm —e— dna? n) B, =o 


[19.10"] 
(nm e+ 4na® n) B= 0 | 

This new system has two solutions: 

(1) Ey = 0; E, => iH, ) n= Ny ) [19.12] 
ni—e—4daain=0, 
ty 
so that (with « sufficiently small) 

— @ 

My met 270 : [19.12] 


This produces a left-hand circularly polarized wave (in ac- 
cordance with Eq. [19.12]). The corresponding index of 
refraction is given by Eq. [19.12’]. 


(2) EH, =0; E, = — iE, , n= Nyy. 


From a calculation analogous to the preceding one, there 
follows 


@ 
Nyy ~ Ve — 20a — . 


Thus, from Eq. [19.1], the rotation of the plane of oscil- 
lation becomes 


2 
X= — Inx (2) L. [19.1'] 


In our theory, conservation of energy follows correctly 
if we make the following assumptions: 


s= 7 Bx(B+ Sak 
4a c 


it 
w= a {eE? + B* — 4aaE-curl E}. 
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The essential feature of this theory is that in the funda- 
mental assumption, Eq. [19.4], we establish a relationship 
between a proper vector and one which has the transfor- 
mation character of a vector product. 

In the case of optically active crystals, the relation be- 
tween D and E must be made more general: 


D; = DnB, . [19.8a] 


Conservation of energy requires that the tensor 7,, be her- 
mitian, 
Qe = Nie ’ 


so that it can be decomposed in the following manner: 


Nix = Ex + Ux » 


with €:,—=€&:, and 0,,——o.;. Consequently, we can write 
further 
On = D firm, [19.8b] 
t 
with 
fes,r Se fit, . 


Equation [19.8a] can also be written with the help of a 
vector product. Let ge be the vector associated with the 
skew-symmetric tensor 9;,. That is, let 


0: = > ecrOn ’ 


k<l 


where ¢é;,, denotes the familiar permutation tensor: 


+1 if i, k, Lis an even permutation of 1, 2, 3; 
Ent={—1 if i, k, lis an odd permutation of 1, 2, 3; 


0 if any two of the indices are equal . 
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Then, Eq. [19.8a] goes over into 
D, = €jn Ly, — i(p x E); 
and the formula corresponding to Eq. [19.8b] reads 
i D fam . 
The equations for isotropic media are re-obtained from these 
two formulas if we set 
fae Oa, Tii—Tem, GH £2 —e0ne 
b. Kuhn’s model [A-7] 


As a model of an optically active molecule we take two 
coupled oscillators which, in what follows, are to be ar- 
ranged in a mirror noninvariant way (see Fig. 19.1). To 


Oscillation of 2 


= 


Oscillation of 1 


Figure 19.1 


simplify the calculation, we assume equal masses for the 
two oscillators. Then, taking the coupling into account, 
the potential energy of the system becomes 


m 
Et — DY (wi ax} Se W393 ap K 213) . 


For free oscillations (in the absence of an electromagnetic 
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field) we obtain 
©, = — wa, — ky, — (1/m)oE,,,/O2, , 
Y2 = — ka, — wey, — (1/m) OF,/OYs . 
With the assumption 
= Aim, y,= Beir, 
there follows 


ee | 
b 


ak + yx(— 03 + 3) = 0 [19.13] 


which, if a nontrivial solution is to exist for 7, and y, 
leads to the determinantal requirement 


(— wo + wi)(— ws + w3) —k? = 0. [19.137] 


The roots w’ and w” will be in the neighborhood of the 
eigenfrequencies @, and w,, corresponding to the uncoupled 
oscillators. Thus, we have the expansions 


k2 
wo! ~ wi — G5 
Ws —s W, 
: [19.14] 
m2 2 ke 
wo" ~ wi + = 
Ws ae WM, 


For the two eigenoscillations of the system, there follows 
from Eq. [19.13] 


= 
Be ato” 19.15 
ye ees si 
On 


If we now assume light incident parallel to the 2 axis, 
then we have a forced oscillation. With the same assump- 
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tion as before, instead of Eq. [19.13] we now obtain 
€ 
(— w*+ w3)a,+ ky, = 7 Ey, 
. , : [19.16] 
kar, + (— w* + 03). = in Ey 


where w is the frequency of the light and all quantities 
are periodic with this frequency. This yields 
[(— * + o4)(— w* + 3) — kJ, 
[(— w*+ @3) F,.— kKEyy] 
. [1926 


S]o 


[(— @*+ wi)(— w?+ ws) — k* ly, 
e 2 2 
os all — w"+ w}) By — kE,s) 


On the basis of Eq. [19.13’], we can write the first square 
bracket in Eq. [19.16’] as 


(w® — w'*)(@* — w"*) . 


Thus, there follows for the electric dipole moment 


- a e?/m 
la 2 ia (w%— w"*) (w*? — w") 
X[(— w* + w3) Ei. — kEy] 
= 2 e?/m » [19.17] 
= 2 Oe (a — wo" 0") 


x [(-— a 3) Ey, oa kE,,| 
P, = > ez, = constant 
k 


If we denote the electric field intensity of the incident 
wave at the center of the interval d by E, then we ean 
expand and obtain 


d ok, 
Be=2,—5 "=... 
By, = B+ Es 


2 oz 
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Thus, 


“ et/m 
P.= wooo w* + 3) E, 
d oH, 2, 2 OH, 
ada 2 Oz See 2 oz 
ae [19 17'] 
a a a aoe 


d ok, dok, 


+ (08+ of)By + (— 08+ ol) = 5 +k 


For anisotropic media, this cannot be further simplified. 
However, for isotropic media we can average over the po- 
sitions of the molecules and can calculate the average 
value of P. The average value of P for all possible orien- 
tations of the oscillators parallel to the x-y plane (denoted 
by P’) is obtained by forming the arithmetic mean of P 
for the original position and P for a position rotated by 90°: 


(- or AE) g 


wi —w}d 0E d 
2 aMae apg eurl E}. 


2 
P= : im nel) 


(w? = w'*)(w? =a) 


The average value (denoted by P) for the possible posi- 
tions of the axes of the oscillators (coinciding with the 
zg axis in the figure) is obtained as an arithmetic mean 
of P’ for the two oppositely directed positions of the axis. 
By calculating P’, making the substitution s>—2, yy, 
g—>—ze and recalculating P', adding the two and taking 
half the sum, we get 


as e?/m 
ia (w* — w'*)(w* — w"*) 
¥ (- gee Oo + “) E+ = curl E [19.18] 


In a completely similar manner we now calculate the 
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magnetic dipole moment, 
€ : 
M= on x REX HK,» 


The magnetic moment must be referred to the same point 
as the electric moment—that is, to the center of the in- 
terval d. There then follows 


é 
M.=—35,5% 
_— Ce ie [19.19] 
MO ye 5 
M,=0 


We now again expand E about the center of the interval d 
and substitute into Eq. [19.19]. (See Eq. [19.17] above.) For 
this, we need go only up to the constant term in the Taylor 
expansion, because Eq. [19.19] is already linear in d and 
we consistently neglect higher-order terms. Thus, there 
follows 


wag OS peal 
open ag 
d e?/m 


== de (w? Sa = w"?) 
x[(- o-oo) #,— Bn] |° [970 


d e?/m 
M,=—-—~— 
4c (w? — w'*)(w* — w") 
x {(— a* + w2) EB —kE,] 
As in the preceding calculations, we must average, first 
over the positions of the oscillators parallel to the w-y plane: 


a e?/m 
M,= 4¢e (on w'*)(w? = ww"? 
—kEH, i— 2 E, 
x[ + $(@; — 3) By] , [19.207] 
Mw -— d e?/m 
( oe 4e (cw? = o'*) (@? ab w"*) 


x[— 3 (w) — wf) B, — kB,] 
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and then over the directions of the axes: 


= 4 e?/m 


= iin orate oy PE [19.21] 


Thus, we see that Kuhn’s model explains the phenom- 
enological quantities «, and «,. By comparison with 
Eq. [19.4], there follows (for isotropic bodies) 


es ee d 
nae (wo* we w'*)(w* = a”) 9 
e?/m d 
a = (w? — o"?) (ep? _ w"*) ka : [19.22] 
2 
Ps e?/m d 


(w? —_ o'*)(@? — w"*) , 4 


To summarize, we want to emphasize that in Kuhn’s 
model the coupling of the oscillators and the position de- 
pendence of the field intensities of the incident wave are 
essential. The decomposition of the effect into a part from 
P and a part from M is not essential; they can be combined 
by introducing the concept of the retarded electric mo- 
ment.* The fact that the finite velocity of propagation of 
the force between the two oscillators has not been con- 
sidered is also unessential. This is true because it can be 
shown that if the forces in the atom are electrical in na- 
ture, this effect is of order of magnitude d? and, conse- 
quently, is to be neglected. Also unessential is the fact 
that we had only two coupled oscillators. For the general 
case, see M. Born, Optik (Supplementary Bibliography). 


20. MAGNETO-OPTICS 


a. Zeeman effect 


In the following, we will study the oscillations of a par- 
ticle in a static magnetic field H. 


See Krooy, Dissertation, Leyden (1936). 
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The fundamental equation is obtained by assuming the 
Lorentz force: 
Se : Ze 
= — , 20.1 
x" + Wo x mo [20.1] 


If we take the z axis parallel to the magnetic field H so 
that H={0, 0, H}, then there follows: 
1. For the 2 component: 


Z2+0j2=0, 
2 —1Cenee. 


This is an oscillation parallel to the magnetic field with 
unchanged frequency @p. 
2. For the x and y components: 


‘s eH . 
B+ wow =— Ys 


i Pr eH F 
Y¥+ooy = ae iE 
With the substitutions 
& =a@ ale ay ’ 


&y = { =ax2—iy, 
we get = 


- ‘ ne. 
& + w= —— ik, 


xs eH .- 
Ey + ooon= te we Ey . 
If we assume 
é, = Ae-‘'t, 


bx — Bees, 
then upon substitution, 
eH 
—w"*+ 05+ —o'=0, 
me 


eH 


2 
— ow" + ow} —— w"=0. 
me 
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This last system of equations is rigorously solvable, but it 
is sufficient to expand the roots w’ and w"’: 


, eH 
ia il 2me [20.2] 
: yeh oad 
Ne as 


The quantity eH/2me is called the Larmor frequency. 

Thus, upon turning on the magnetic field, the oscillation 

of frequency , has been split into 
1. an oscillation parallel to H (the so-called x oscillation), 
with unchanged frequency a; 
2. two oscillations, one right-hand and one left-hand cir- 
cularly polarized, both perpendicular to H (the so-called 
o oscillations), whose frequencies are symmetrically dis- 
placed from w, by an amount eH/2mc. The oscillation with 
the higher frequency (for ¢ positive) is left-hand circularly 
polarized [A-2]. 

For the case of electrons emitting light, this splitting of 
the oscillations gives rise to a-splitting of the spectral lines. 
If observations are made transverse to the magnetic field, 
then the z oscillation appears linearly polarized parallel 
to H and the o oscillations appear linearly polarized per- 
pendicular to H. (By the polarization direction, we mean 
here the direction of the electric vector.) If observations 
are made along the direction of the magnetic field, then 
the z oscillations are not seen at all, while the o oscilla- 
tions appear circularly polarized with the higher-frequency 
oscillation being left-hand circularly polarized. From this, 
however, it follows that the radiating charge must have a 
negative sign because, otherwise, if the observation direc- 
tion were reversed, the higher-frequency oscillation would 
be associated with the left-hand circular polarization [A-2]. 

The Zeeman effect permits the determination of e/m (in- 
cluding the sign); the value thus obtained is that for 


electrons. 
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The simple Zeeman splitting considered above is not the 
most general case. For many spectral lines more compli- 
cated splittings with several a and o oscillations appear; 
this is the so-called anomalous Zeeman effect. This effect 
is due to the spin of the electron. However, the correspon- 
ding theory can be treated only with wave mechanics. 


b. Faraday effect 


This effect consists of the rotation of the plane of oscil- 
lation of light in a magnetic field H. The refractive index 
is different for right-hand and left-hand circularly pola- 
rized light. 

1. Light rays parallel to H. The fundamental equation 
becomes 


#+ tx = — &xH+—E. [20.3] 


Now, 
P= Nex, 


where N is the number of particles per unit volume. We 
choose H parallel to the z axis. Let the incident wave have 
angular frequeney w, so that x, which will execute forced 
oscillations, also has this frequency. We make the follow- 
ing substitutions: 


E,+ ik, = Eyn ’ 
FASE 1D, Sar pa, 


from which follows 
eH e2 
(w" +0840) Pin = NT fen : 


Upon neglecting the corrections between E’ and E we 
have Eq. [16.7]: 


4nP = (n®?—1)E. 
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Thus, 


a 4at(e2/m).N 
wh — w+ (eH/me)w 
= 4(e2?/m)N 


ws, — w* — (eH/me)w’ 
so that 
1 
wo — w* + (eH /me) w 


2 2 e 


il 
. {= aiimara ; 


With the assumption that w does not lie close to the eigen- 
frequency w,, this can be simplified (mean value theorem): 


e? il eH 
m (Wyg—a")" me [20.4] 
e? al w eH , 


Se ae gy | een Aaa 
ie ae (Oey a me 
Now, however, n, and n,, are by hypothesis just the re- 
fractive indices for the two circular components of the 
linearly polarized incident wave. Thus, in accordance with 
Eq. [19.1], the angle of rotation y of the plane of oscil- 
lation becomes 
y=vlH, 

with 

__ 4n(e?/m)N e w? 

SS : 20.5 

~ (w3— w*)? me? 2n ad 
The quantity v is called Verdet’s constant. For the case 
of several particles in a molecule one must, of course, use 


the following sum for ?: 


- o= 5 ele eee [20.6] 

)? mc? 2n° 
2. Transverse incidence. In this case birefringence will be 
observed: one oscillation parallel to H and one perpendi- 


cular to H. 
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From the technical point of view, the Faraday effect is 
much easier to detect than the Zeeman effect, although 
from the theoretical standpoint the Zeeman effect is sim- 
pler. The Zeeman effect consists merely of free oscillations 
of electrons which are influenced by a static magnetic field, 
while for the Faraday effect forced oscillations under other- 
wise identical conditions occur. 
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Appendix. Comments by the Editor 


[A-1] (p. 23). This is a Legendre transformation from 
x, x° to p, p® with dS = p-dx—p°-dx°®. The latter form 
is a special case of Eq. [3.12]. 


[A-2] (pp. 31, 68, 145). The case a= 0 gives a clockwise 
rotation which by convention is called right-hand circular 
polarization; b = 0 gives a counterclogkwise rotation called 
left-hand circular polarization. Convention fixes the sense 
of rotation as seen when looking against the direction of 
light propagation. It also fixes the plane of polarization 
to be spanned by the magnetic field vector and the wave 
vector. See M. Born, Optik, p. 24 (Supplementary Bib- 
liography). 

[A-3] (p. 34). With the advent of high-power lasers 
(coherent light sources) nonlinear optics has become a 
fascinating reality (excitation of harmonics, light-light 
scattering, etc.). See, for example, N. Bloembergen, Non- 
linear Optics (Benjamin, New York, 1965). 


[A-4] (p. 62). It follows from Eqs. [7.8], [7.8'] with a, = f, 
= P,=0, « = sing that 


pa 

m1 pik(RoptRip) 

u = constant 3S (f erinbeinede , 
p=0 RopRrs 


where a is the height of the steps and R,, = R,-+ npd, 
Ry, = R, — pdcosg, so that, for small oy, Ry, + Ry = Ri 
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+R, + pd(n—i1). (The denominator may be approxi- 
mated by R,R,.) Integration and summation yield 


1 — etkasing | _ pikmid(n—1)—asin 9) 
uw = constant 


sing 1 — giklain—D-asing) * 

The first factor vanishes for (a/A) sing = +1 which for 
a> A= 2ax/k limits the refraction angle g to small values. 
The principal maxima are given by the simultaneous 
vanishing of numerator and denominator of the second 
factor, 


m[d(n —1)—asing]=N,/, N, integer ; 


d(n—1)—asing =WNA, N integer. 


Within the interval | (a/A)sing|<1 there is only one 
principal maximum N= (d/A)(n—1). The zero next to 
it is displaced by Aq as given by 


tg a 


1 


On the other hand a change of wavelength corresponds to 
a displacement 3p given by 


dn d 


The resolving power is determined by the condition 
Ag = 8 and is given by 


[A-5] (p. 108). The situation is as presented in fig- 
ure A.l. See M. Born, Optik, p. 238 (Supplementary Bib- 
liography). 
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[A-6) (p. 113). This is the f-sum rule of Thomas and Kuhn 
(see M. Born, Optik, pp. 468, 525, Supplementary Bib- 
liography). In quantum mechanics the f-sum rule is a 
consequence of the commutation relation between the 


m,-cone 


Figure A.1 


coordinates «, and the velocity components 2%,, 2, %, 
ee ao Oy: 

[A-7] (p. 138). W. Kuhn, Z. physik. Chemie B4, 14 (1929). 
See M. Born, Optik, p. 527 (Supplementary Bibliography). 
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Cornu spiral, 64 

Crystal optics, 16, 84-108 
Crystals, 20 
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Current density, 109 
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Damping constant, 79 
Damping force, 117 
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Fraunhofer, 53 
Fresnel, 53, 62-66 
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by slit, 54-55 
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Dispersion law, 34, 35, 38, 117 
for x rays, 128 
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Dual equations, 95, 96 

Dual variables, 92 
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Eikonal equation, 8 
Electrons, 113, 119, 127-131, 145 
Electron, spin of, 146 
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Energy conservation, 16 
for electromagnetic waves, 69, 85, 
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Energy density, 
of electromagnetic field, 85 
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total, 116 
Energy transport, direction of, 93 
Euler, L., 28 
Euler’s homogenity relation, 16, 99 
Extremal curve, 5 
Extremal principle in canonical 
form, 15 
Extremal property, 2. 


Faraday effect, 146-148 
Fermat’s principle, 3, 4, 6, 7, 14 
Fluctuations, 

density, 124 

theory of, 126 
Focal plane, 61 
Fourier inversion formula, 35 
Fourier’s theorem, 31, 35 
Fraunhofer diffraction, 53, 54-62 
Fresnel, A. J., 28 
Fresnel diffraction, 53, 62-66 
Fresnel formulas, 6 
Fresnel integral, 63 
Fresnel mirror-experiment, 43 
Fresnel’s formulae, 

for absorbing media, 79-83 

for nonabsorbing media, 73-77 
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Larmor, 145 
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Grating, 57, 58, 60, 61 
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Green’s theorem, 46 
Group velocity, 37, 98, 115, 116 
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Hamilton’s theory, 
homogeneous form of, 14 
inhomogeneous form of, 21 
Harmonic oscillation, 28 
Hertz oscillator, 123 
energy radiation of, 124 
Hertz vector, 122 
Huygens, C., 106 
Huygens’s principle, 9, 20 


illuminance, 26, 42 
Image plane, 61 
Image space, 1] 
Index of refraction (refractive 
index), 2, 40, 111 
of crystals (double root), 88, 
100, 101 
in Faraday effect, 146 
of metals (complex), 77-83, 120 
in optical activity, 133, 136 
variable, 4, 28 
Intensity, 29, 43, 57, 58, 74, 125 
incident, 75, 130, 131 
of illumination, 26 
reflected, 75, 80 
refracted, 75 
of scattered x rays, 130 
specific, of light, 26 
time-averaged, 30 
Intensity maxima, 44 
Intensity ratio, 65 
Interference, 28, 42, 44, 45, 67 


Kepler, law of, 25. See also 
Lambert, law of 
Kirchhoff-Clausius law, 12 
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Kirchhoff, solution of, 45—53 

Kuhn, f-sum rule of Thomas and, 
152-153 

Kuhn’s model, 138-143 


Lagrangian multipliers, method 
of, 3, 15, 90 

Lambert, law of, 25, 26 

Larmor frequency, 145 

Lasers, 151 

Legendre transformation, 151 

Lenses, 13, 53 

Light, 
circularly polarized, 133, 134, 146 
corpuscular theory of, 28 
elliptically polarized, 77, 81 
linearly polarized, 75, 76, 77, 80 
polarized, 67 
unpolarized, 67, 75, 131 

Light intensity, specific, 26 

Light vector, 67, 68 

Lorentz force, 144 

Lorentz’s theory of scattering, 124 


Magnetic field (static), 143, 146 
Malus’s law, 9 
Mass absorption coefficient, 128 
Maxwell equations, 68, 110, 134 
in crystal optics, 84 
for periodic processes, 71, 77, 
86, 121 
Mechanical momentum, 16 
Mechanical problem, associated, 14 
Mechanical waves, 34 
Media, 
absorbing (metals), 77-83, 119-121 
anisotropic, 141 
dispersion-free, 39 
homogeneous, 1, 124 
inhomogeneous, 41 
isotropic, 1, 34, 39, 41, 42, 45, 71, 
116, 138, 141 
nonabsorbing, 71-77 
Metals. See Media, absorbing 
(metals) 
Michelson’s transmission echelon. 
See Grating, step- 
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Microscope, 59 
Mirror, 
Fresnel, 43 
perfect, 82 
Molecules, not mirror-symmetric, 
133 


Newton’s rings, 45 

Nodes of E and H in standing 
waves, 83 

Nonlinear optics, 151 


Object, 

illuminated, 60 

luminous, 59 
Object space, 11 
Optical activity, 133 

Optical axes, 90, 101 
Optical constant for metals, 82, 120 
Optical instruments, 62 
Optically active crystals, 137 
Optically active molecule, 138 
Optically uniaxial crystals, 90 
Optical path, 3, 4, 6, 107 
Oscillation, 

damped, 118 

harmonic, 28 

a (parallel to H), 145 

o (perpendicular to H), 145 
vectorial, 30 
Oscillator, 

Hertz, 123 

quasi-elastic, 110 
Oscillator strength, 113 


Parseval’s formula, 33, 36. See also 
Completeness relation 
Pedal surface, 97 
Phase, 
continuity of, 40, 72 
propagation of, 93 
Phase velocity, 34, 40, 68, 71, 86, 
88, 110 
Photometric brightness. See Light 
intensity, specific 
Plane, 
of incidence, 2, 72, 73, 79, 80, 81, 82 
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Plane (continued) 

of polarization, 68, 76, 81, 83, 151 
Polarization, 67, 77, 79 

circular, 133, 145, 151 

electric, 111, 124 

plane of, 68, 76 

for plane waves, 69, 79 
Potentials, 

scalar, 121 

vector, 121 
Poynting vector, 74, 77, 84, 91, 

93, 114 

Principal angle of incidence, 81 
Principal axes, 85 
Principal dielectric constants, 86 
Principal maxima, 58, 59, 61, 152 
Principal section, 88, 89 


Quantum-mechanical model, 120, 
131 

Quantum (wave) mechanics, 146, 
153 


Radiation, energy, 124 
Radiation condition, 47 
Radiation damping, 121 
Radiator, isotropic, 26 
Ray, 
extraordinary, 89 
ordinary, 88, 89 
parallel, 53 
reflected, 2, 75 
refracted, 2, 75 
Ray axes, 105 
Ray cone, 104 
Ray direction, 20, 93, 97-107 
Ray equation, 10, 16, 19, 22 
Ray index, 91, 97, 98 
Ray surface, 98 
Ray variables, 91 
Ray vector, 91 
Ray velocity, 91 
Reflection, 
law of, 2, 4, 41 
total, 75, 76, 77 
Refraction, 
double, 108 
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index of (refractive index), 2, 4, 
28, 40, 79, 86, 111, 120, 133, 136, 
146 

law of, 2, 4, 41 

Refractivity, molar, 113 
Resolving power, 59-62, 152 
Rotation of plane of oscillation, 

angle of, 133, 147 

in Faraday effect, 146 

in optical activity, 133, 136 


Scattered energy, 125 
Scattering, light-light, 151 
Scattering angle, 132 
Scattering coefficient, 126, 127 
Screens, 46, 48, 49 
Shadow boundary, 51, 65 
Sine condition, 12, 60 
Singular directions, 99, 105 
Shit, 54, 63 
Snell’s law, 1, 7 
Spectral lines, splitting of, 145, 146 
Statistical average, 132 
Statistical mechanics, 124, 126 
Structure factor, 130 

in quantum theory, 132 
Superposition principle, linear, 

29, 34, 35, 38 
Surface, 

of constant phase, 86, 107 

pedal, 97, 98 

ray, 98 

wave, 20, 97 
Surface wave, 77 


Tensor, 

permutation, 137 

skew-symmetric, 134, 137 
Thomas and Kuhn, f-sum rule of, 
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Thomson, J. J., 128 
Total reflection, 75, 76, 77 

limiting angle of, 76 
Transversality condition, 18 
Transversality of light waves, 28, 67 


Variational principle, 11, 16, 17, 21 
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Vectors, 

axial, 134 

polar, 134 
Verdet’s constant, 147 


Wave equation, 38 
general solution of, 39 
Kirchhoff’s solution of, 48, 49 
time-independent, 71 
Wave, 
damped, 79 
electromagnetic, 69, 74 
harmonic, 43, 45 
incident, 71, 72, 73, 74 
mechanical, 34 
outgoing, 47 
plane, 33, 39, 69, 71, 72, 80, 86, 
107, 135 
reflected, 72, 73, 74, 76, 82 
refracted, 73, 74, 76, 79, 80 
scalar, 67 
spherical, 39 
standing, 82 
surface, 77 
transverse vector, 67 
Wave concept, 28 
Wave cone, 105 
Wave kinematics, 28, 38 
Wavelength, 42, 45, 59, 124, 152 
Wave normal, 20, 34 
Wave packet, 35, 37, 45, 93 
Wave surface, 20, 97 


X rays, 111, 124, 128 
scattering of, 127-133 


Young, T., 28 


Zeeman effect, 143-146 
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